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HAUPTAUFSÄTZE 


Zur Theorie und Konstruktion der Geschoßflugbahn. 
Von @. Schweikert in Charlottenburg. 
Es wird ein neues, einfaches und doch strenges und allgemeines Verfahren 
zur graphischen Konstruktion der Geschoßflugbahn entwickelt. Unter Zu- 
grundelegung einer zusätzlichen, normalen Widerstandskomponente werden 
strenge und allgemeine ballistische Hauptgleichungen abgeleitet und Ver- 
fahren zu ihrer rechnerischen und mechanischen Integration mittels einer 
neuen Integriermaschine sowie zur experimentellen Ermittlung der Wider- 
standsgesetze angegeben. 
Einleitung. 

Wir setzen bei unseren folgenden Untersuchungen zunächst voraus, wie dieses in der 
Theorie der äußeren Ballistik üblich ist, daß die Resultierende des Luftwiderstandes des Ge- 
schosses in die Tangentenrichtung der Geschoßbahn fällt, nnd daß somit senkrecht zu der 
Tangente, d. i. in Richtung der Normale der Geschoßbahn keine Komponente auftritt. Ob- 
schon dieses streng nur dann gilt, wie man leicht durch eine einfache Überlegung zeigen 
kann, wenn die Symmetrieachse der Geschoßform, also praktisch meistens die Geschoßachse 
selbst durch den Schwerpunkt des Geschosses geht und genau in die Tangentenrichtung der 
Geschoßbahn fällt, was ja auch bei den durch Drall stabilisierten Geschossen bekanntlich nie 
zutrifft, hat man doch bisher in der ‚Literatur die genannte einschränkende Bedingung nicht 
genauer festgelegt und formuliert. 

Wenn die Symmetrieachse der Geschoßform mit der Bahntangente zusammenfällt, so 
gibt es zu jedem Flächenelement der Geschoßoberfläche ein genau symmetrisch gelegenes 
Flächenelement, das den genau gleichen Luftwiderstandskräften nach Größe und Richtung 
unterworfen ist. Zerlegt man nun den an jedem Flächenelement angreifenden Luftwider- 
standsvektor in zwei Komponenten, von denen die eine in die Tangentenrichtung, d. i. also 
in unserem Falle auch in die Richtung der Symmetrieachse fällt, während die andere Kom- 
poneni. die dazu senkrechte Richtung der Normale der Geschoßbahn hat, so heben sich eben 
deshalb, weil die Bahntangente nach Voraussetzung auch Symmetrieachse der Geschoßform 
ist, die normalen Komponenten in den einzelnen Flächenelementen gegenseitig auf, so daß 
nur die tangentialen Komponenten übrigbleiben. Nur in diesem besonderen, einfachen Falle 
hat dann die in der Ballistik allgemein gebräuchliche Zerlegung des Luftwiderstandsvektors 
(W) in die horizontale Komponente (W,) und die vertikale Komponente (W,) in der ein- 


fachen Form: W,=W-cosd; W,=W-sind 


Gültigkeit, wobei # der Neigungswinkel der Bahntangente gegen die horizontale x-Achse ist. 
Ferner gilt auch nur in diesem Spezialfalle die bekannte Gleichung für die zentrifugale 
Beschleunigung des Geschosses, wie sie allgemein in der Literatur zu finden ist: 
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v’:0=—g:-cosd, 


wobei v die Geschoßgeschwindigkeit, o der Krümmungsradius der Geschoßbahn und 9 die 
Erdbeschleunigung ist. 

Allgemein muß man demgegenüber beachten, daß je nach der Geschoßform und der 
Lage der Symmetrieachse dieser Form zur Bahntangente nicht nur eine tangentiale Kompo- 
nente Wr, sondern auch eine normale Komponente Wy des Luftwiderstandsvektors auf- 
treten kann’). 

In Bild 1 bezeichnet & die Richtung der horizontalen 
Abszisse, T die Richtung der Bahntangente, N die Richtung 
der Bahnnormale, S die Richtung der Symmetrieachse des 
Geschosses und @ die Richtung der Erdbeschleunigung in 
dem betrachteten Bahnelement. Der Winkel zwischen x 
und 7 sei wie üblich mit # bezeichnet, derjenige zwischen 
T und $ hingegen mit a. Dieser letztere ist also der so- 
genannte Anstellwinkel der Geschoßachse oder genauer der 
Symmetrieachse der Geschoßform gegen die Bahntangente. 

Je nach der besonderen Form des Geschosses und der 
Größe des Winkels «a bekommt man stets außer der Luft- 
widerstandskomponente ‚Wr in Richtung der Bahntangente noch eine Widerstandskomponente 
W,„ in der dazu senkrechten Richtung der Bahnnormale. Letztere kann nun je nach der be- 
sonderen Lage des Anstellwinkels und der besonderen Geschoßform entweder nach oberhalb 
der Tangente gerichtet, d. i. positiv sein und einen Auftrieb des Geschosses bedingen, oder 
aber unterhalb der Tangente liegen, d. i. negativ sein und einen Abtrieb des Geschosses 
zur Folge haben... Die Möglichkeit der seitlichen Abweichung aus der Mündungsebene der 
Geschoßflugbahn lassen wir dabei außer acht. 

Die beiden Komponenten Wr und W, lassen sich natürlich wiederum zu einer einzigen 
Resultierenden des Luftwiderstandes zusammensetzen, die mit der Bahntangente den Winkel y 
bilden möge. Dann gelten die, einfachen Beziehungen: 


ty =p,: y=arctg Wr 


Für die horizontale bzw. vertikale Komponente des Luftwiderstandsvektors W, und W, aber 
gelten allgemein die Gleichungen: 


W;= Wr: cos® + Wy-sin®, 
W,= Wr: sind — Wy:cos®. 


Bezeichnet man schließlich den Winkel: 








+y=d, 


und den Betrag des resultierenden Widerstandsvektors mit W,, so läßt sich statt dessen auch 
einfacher schreiben: 

W,=W,:cos®%,, 

W, —W, ” sin d, 
in offenkundiger Analogie zu den bisher üblichen Formulierungen. Dabei ist aber der Winkel 
9%,=d+y an die Stelle des Winkels # getreten, wobei y der Winkel zwischen der Luft- 
widerstandsresultierenden und der Bahntangente ist, und nur mit dieser Festsetzung gelten 
die vorstehenden einfachen Beziehungen allgemein. Auch dieser Umstand ist bisher in’ der 
ballistischen Literatur nicht genauer formuliert worden. 


Schließlich tritt noch für die zentrifugale Beschleunigung des Geschosses an die Stelle 
der gebräuchlichen Gleichung: 


v:0=—g:cosd 
allgemein die genauere Beziehung: 
"_w ä (5) 
TEEN 
au. . r u u 
wie sich. aus Bild 1 ohne weiteres ablesen läßt, wobei di die Winkelgeschwindigkeit der 
Bahntangente ist. 


1) Vgl. C. Cranz: Lehrbuch der Ballistik, 1. Bd. Berlin 1935, $ 17, S. 108, sowie Ergänzungshand Berlin 1936, 
% 12, 8.32; Th..Vahlen: Ballistik, Berlin 1942, S, 16 ff. 
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Wir beschränken uns im folgenden zunächst auf den einfacheren Spezialfall, daß die 
Symmetrieachse der Geschoßform mit der Tangente der Geschoßbahn zusammenfällt und somit: 
Wı =0; y=P0; 9, =d 
ist und formulieren anschließend die genaueren Gleichungen, welche auch der Normalkompo- 

nente des Luftwiderstandes Rechnung tragen. 


8 1. Graphische Konstruktion der Geschoßbahn. 
Wir gehen aus von den beiden bekannten Differentialgleichungen der Geschoßbahn: 
v’ a do\? 
r = —9-cosd=90: (2) 
für die normale Richtung und: 
dv 
dt 
für die tangentiale Richtung, wobei die letztere Gleichung auch allgemein für beliebige An- 
stellwinkel der Symmetrieachse zur Bahntangente Gültigkeit hat, da ja der Beschleunigungs- 
vektor immer die Richtung der Bahntangente hat. 
Aus der ersten der vorstehenden Gleichungen erhält man sofort: 


— Wr —g- sind 


d9 dd v —g:-cosd g 
m29.+,!% = — _— ı[ — 
vier we v u.’ 
wenn wir noch: ö 
u 
cosd r 


setzen. Die Größe v, kann man leicht graphisch gewinnen, indem man auf dem Geschwin- 
digkeitsvektor in seinem Endpunkte eine Senkrechte errichtet, die auf der Abszissenachse x 
den Wert v, abschneidet. Verbinden wir dessen Endpunkt ‚mit dem Endpunkt des Erd- 
beschleunigungsvektors g und bezeichnen den Winkel an der Spitze mit ß, so gilt: 
g me v| 
tg ß Ur | Q | 

Man gewinnt also den Krümmungsradius in dem betrachteten Bahnpunkte sofort, wenn man 
in Abstande v eine Parallele zur x- Achse zieht und von deren Schnittpunkte mit der zuletzt 
gezogenen Verbindungslinie, d. i. dem freien Schenkel des Winkels 3 eine Senkrechte auf 
die Abszissenachse fällt. Der Abstand des Fußpunktes derselben von dem Scheitelpunkte 
des Winkels £ gibt die Größe des Krümmungsradius o (vgl. Bild 2). 

Wir errich- 7 
ten sodann in N u h 
dem betrachteten D, 9 „fh 
Bahnpunkt Peine 2 
Senkrechte auf p # 2 
dem Geschwin- 
digkeitsvektor v, y ® fi 
d.h. auf der Bahn- 
tangente und tra- v 
gen auf dieser o ER 
ab, so erhalten e 
wir den Krüm- 
mungsmittelpunkt Bild 2 
O0 der Geschoß- 
bahn im Punkte P. Schlagen wir um O einen Kreis mit o, so 
liefert dieser den Krümmungs- oder Oskulationskreis,: der den 
Verlauf der Flugbahn im Punkte P für den Bereich des Winkel- 
elementes d# darstellt. Wir ziehen also einen zweiten Radius 
von O aus, der den kleinen Winkel 4% mit O P bildet und uns 
den neuen Flugbahnpunkt P, liefert. 

Die Senkrechte in P, auf OP, gibt die Richtung des neuen 
Geschwindigkeitsvektors v,. Um auch seinen absoluten Betrag zu 
ermitteln, benutzen wir die zweite der obigen Grundgleichungen: 


dv z ; 
de Wr—g:-sind. 


N 































Dabei setzen wir voraus, daß der Geschoßwiderstand: 
Wr=e-ö(y):v’-K(v, y) 
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durch eine Tabelle oder Kurve als Funktion der Geschwindigkeit v und der Höhe y für be- 
liebige Luftzustände (Temperatur, Wind und Feuchtigkeit) gegeben ist. Desgleichen kann 
die Erdbeschleunigung g tabellarisch als Funktion der Höhe y festgelegt sein. Die Strecke 
g-sin® gewinnt man leicht an Hand von Bild 2, indem man vom Endpunkt G@ des Be- 
schleunigungsvektors g das Lot auf den Krümmungsradius O0 P fällt. Die Länge dieses Lotes 
ist g-sind. Wir formen nun die obige Gleichung um, indem wir beiderseits durch ® divi- 


dieren und erhalten so: 
dv  _Wr+g:sind 


v-dt v 


Die Strecke »-dt ist aber nichts anderes wie die Länge des Bogens unseres Krümmungs: 
kreises zwischen P und P,. Wir gewinnen somit dv leicht durch eine Ähnlichkeitskonstruktion. 
Wir tragen die Strecke v:dt auf dem Geschwindigkeitsvektor v von P aus ab, oder ersetzen 
einfacher die Strecke v-dt durch das zugehörige Tangentenstück PP,, das der Radius OP, 


auf der Tangente (v) in P, abschneidet. 


Tragen wir also die Strecken Wr und g-sin®ö nacheinander auf der Abszissenachse 
von P aus ab, verbinden den Endpunkt mit dem Ende des Geschwindigkeitsvektors » und 
ziehen durch den Endpunkt P, von v-dt eine Parallele zu dieser Verbindungslinie, so 
schneidet diese auf der Abszissenachse die Strecke Av ab. Die Größe des neuen Ge- 
schwindigkeitsvektors v, im Punkte P, ist dann sofort bestimmt: 





v=v— dev. 


Wir tragen somit auf der Senkrechten, die in P, auf‘OP, errichtet ist, die Strecke 
v— Av ab und haben damit den neuen Geschwindigkeitsvektor in P, erhalten. Nunmehr 
kann von diesem aus die gleiche Konstruktion, wie vorstehend von v ausgehend beschrieben 
wurde, wiederholt werden, um einen neuen Flugbahnpunkt P, und den zugehörigen Ge- 
schwindigkeitsvektor v, zu erhalten?).. Fällen wir von P, auf die «-Achse das Lot und 
bezeichnen seinen Fußpunkt mit F, so ist PF=4Ax und P,F=4y. 


In Bild 2 bestehen also folgende Beziehungen: 
P&6=g; PT=v; PiA=v; PC=DE=v,; AD=PO=o0 XPAG=Bß; 
PP,=v-4t; APOP,=49; PB=g:snd+Wr; PH=4v; 
PT,=u,=v—-Av; PF=4za; P,F=4y. 


Man kann diese Methode noch leicht in dem einleitend dargestellten Sinne verallgemeinern 
und der Normalkomponente Wy des Luftwiderstandes Rechnung tragen, wenn man in der 
obigen Definitionsgleichung für den Winkel $ den Wert g der Erdbeschleunigung durch den 


Ausdruck: 





Wy 
Te0sd 9 
d. h. durch einen reduzierten Wert g, ersetzt: 
I _|ev| 
wi=,-|,| 


Ferner kann man für eine genauere Konstruktion von Fernbahnen der Veränderlichkeit 
des Vektors g nach Richtung und Betrag mit der Höhe und der Zentrifugalkraft bzw. Coriolis- 
kraft in bekannter Weise Rechnung tragen, indem man den entsprechenden Vektor @ der 
Erdbeschleunigung in den einzelnen Bahnpunkten (P,) der Geschoßbahn in seiner Neigung 
gegen die senkrechte Ordinatenachse veränderlich zeichnet entsprechend der jeweiligen 
Änderung der geographischen Breite @ und des Korrekturwinkels: 


n’: R-sin’o, _n?’-R-sin’p 
29 29 


wobei R der Erdradius, » die Winkelgeschwindigkeit der Erddrehung, g, die Erdbeschleunigung 
im Pol und , der Winkel zwischen der Äquatorebene und dem Erdradius zum Anfangs- 


punkt (P,) der Geschoßbahn ist: 





sind= 


P=Mm+tB. 





?2) Die hier beschriebene neue Methode ist grundsätzlich verschieden von den bekannten graphischen Verfahren 
von C. Cranz und R. Rothe sowie von Th. Vahlen, hingegen bestehen einige Berührungspunkte mit den ältesten 
graphischen Methoden von Poncelet: Lecons de m6canique industrielle. Metz 1828/29, S. 55, und Didion: Trait6 de 
ballistique, Paris 1848, S. 196. 
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Es bedarf kaum der Erwähnung, daß man die hier beschriebene Methode zur stück- 
weisen graphischen Konstruktion der Geschoßbahn auch ohne weiteres mit Hilfe der ab- 
geleiteten Gleichungen zur rechnerischen Festlegung der Geschoßbahn durch einzelne Tabellen- 
werte verwenden kann, und zwar erhält man zunächst”die Geschoßgeschwindigkeit v als 
Funktion des Neigungswinkels 9, während man die Zeit t sofort aus der Beziehung: 

NE 


® 
gewinnt. Schließlich bekommt man die «- und y-Ordinaten aus den Gleichungen: 


ze cosß; Axz=v-At-cosd=g-Ad-cosd; 
I. sind; Ay=v-4t-sind=o-A%-sind. 


8 2. Ableitung einer neuen Hauptgleichung. 


Wir wollen nunmehr eine neue, wichtige und interessante Form der Hodographen- 
gleichung ableiten. Dazu g=hen wir von den beiden Komponentengleichungen der Kraft in 
tangentialer Richtung, d. i. in Richtung des Geschwindigkeitsvektors (v) im Schwerpunkt des 
Geschosses und in normaler Richtung, d. i. in Richtung des zentrifugalen Beschleunigungs- 
vektors (v?:o) aus. Für die tanggntiale Beschleunigung erhält man sofort die Beziehung: 

dv 

dr rg sind En re u a TE A 
während sich für die normale Beschleunigung auf Grund der. bekannten Beziehung für die 
Zentrifugalkraft mit dem Krümmungsradius o der Geschoßbahn und der Winkelgeschwindigkeit 
d#:dt die Gleichung ergibt: 

v? dd 

EEE er ae SE ee a (2). 
Darin ist Wr die tangentiale Komponente und Wy„ die normale Komponente des Luftwider: 
standes, erstere in Richtung des Geschwindigkeitsvektors, letztere in der der Erdbeschleunigung 


entgegengesetzten Richtung positiv gerechnet. 
Aus den vorstehenden Gl. (1) und (2) erhält man durch einfache Umformung die Be- 


ziehungen: 


g-sind=—(Wr-+v); vg. Ba le En St ee 
g:cosd=Wy—v-W; 0-5. en ee ee EEE 


Als weitere dritte Gleichung tritt hinzu die Beziehung für die vertikale Geschwindigkeits- 
komponente: 





dy } v 
=; =v:-ind =——.(W: Biere. 
dt t g (Wr+ev) (3) 
Die vorstehende Gl. (1) differentieren wir noch einmal nach t und erhalten: 
dv dWr d ' dd 
= a nd | 
N (4). 
dWr LE dv; 1. De yv-sind+W,:v' | 


Wir erweitern die Gl. (4) mit v und die Gl. (2) mit g-cos®# und subtrahieren beide von- 
einander, so ergibt sich: 


Wy : 9: c0os# — 9’: co? d+v-(v""+Wr+g'-sind)=0; 
‚_dg_dg dy en 
g- 1 = r >= v-sind; „=-7% 
Hierin führen wir für Wr und g’ die oben abgeleiteten Beziehungen ein: 
en :v- sind — W,- vv” —v:v"=0 
” . (6). 
(+9, :%)- co®#— Wy-g-cosd + W,-— r he +W)-W,-0:v 0-0" 00 
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Dieses ist eine quadratische Gleichung für cos®#, die wir nach g:cos®# auflösen: 




















gs — Wr Zi 
149, 6) 5 
ir Er a MEINES er 
/ 1 u Knien vv’ +gy U? — ze "+ Wr) 
Va el 
9 9 


Da der Winkel d# nur zwischen +7 und -5 variieren kann, so ist g-cos» stets positiv 


und mithin kann nur das positive tie vor dem Wurzelausdruck praktische Bedeutung 
haben: 








1 Wy v ch , Br ERBE 
Ro Ze en 4 at Ze +W,-v Wr -W,-Wr +3.) (8), 
De er Be Er ie ae are Ki 
dabei ist: 
7 e pi, De ZB F BE „m , 
&E=Qy' "ig | 7a) 1 en (10). 


Wir quadrieren beide Gleichungen und addieren sie: 


r 1 Wwx 7 2) (7 ‚Wr z 
—( rtv‘) . Tree age rn Ir" +{m „-W, 'g v”+ Iy- Ww, 2) ® 


Wr yIW m rm). 
oe gti &g o er "+9 FRE 
Zur Abkürzung der Schreibweise fassen wir vorübergehend die Glieder, welche die Ab- 


leitungen der Widerstandsfunktion Wr enthalten, durch eine Bezeichnung zusammen, indem 
wir setzen: 

















(m. - #0). +9 Wr)-0=W, ER a A 
9 I 
Wir schaffen den Wurzelausdruck auf eine Seite und quadrieren die ganze Gleichung: 


1 Wwx [23 
Turptir, -w"+ Wo) erh 


Vorstehender Ausdruck läßt sich wiederum in folgender übersichtlicher Weise zusammenfassen: 





[9’-(Wr+vW)P—2- ))-[a’-(Wr+v')]+ w”+W,)’=0. 








r2 


2 PB. ne 9 2 T, na 6 
Wr Wr 02 


Damit haben wir eine Differentialgleichung zweiter Ordnung und zweiten Grades für 
die Geschoßgeschwindigkeit v als Funktion der Zeit t gewonnen. Die darin enthaltenen 
partiellen Differentialquotienten des Widerstandes W, und der Beschleunigung 9, nach der 
Höhe y können nach bekannten Gesetzmäßigkeiten mit ausreichender Genauigkeit als konstant 
angesehen werden. Dabei ist, wie schon oben erwähnt: 


„= —2: = 308: 10-* sec-?; 


R=Erdradius, W,=— 1,1: 10°*. W7,, wobei Wr, der Widerstand bei der Luftdichte ö, ist. 
Die Normalkomponente W, wie auch die Tangentialkomponente Wr des Luftwiderstandes 
sind allgemein noch Funktionen des Anstellwinkels « der Symmetrieachse der Geschoßform 
gegen die Bahntangente. Wr ist überdies als Funktion der Geschwindigkeit v als gegeben 
vorausgesetzt. 

Da die vorstehende Differentialgleichung die unabhängige Variable t nicht explizit ent- 
hält, so läßt sich ihre Ordnung erniedrigen durch die Substitution: 


gr ‚de, d r 
di=—r Ve Te SEa3. ra 











wobei der Querstrich - Ableitung nach v bezeichnet. Damit ergibt sich folgende Differential- 
gleichung erster Oränung: 
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:_(W, RER ZESP DE IRB. DER an, u 2 _4v 
RR Mr -7z, erde Veeet; 9e (14); 
worin: Eu pie 

u u. W, W x 

W. z= (m, —_ Ey . 0) a Sun (77 Wr 5} ® 
ist. Durch einfache Umformung gewinnt man schließlich folgende Form. der Differential- 
gleichung: 





dp @-(Wr+p?” Wyx — 
Hr VPE FR 
ie ie (15), 
5), +7) 
-(.9-5)2+9,—-—.,)=0 
7 79) u" 
dy_ v (Wr 
1--4-(7 +1 ee a 











Diese beiden simultanen Differentialgleichungen erster Ordnung gelten vollkommen streng 
ohne jede Vernachlässigung für beliebige Flugbahnen des Geschosses. Man kann die Auf- 
lösung derselben wesentlich vereinfachen, wenn man die zusätzliche Abhängigkeit der Größen 
Wr und g von y bei dem Integrationsprozeß vernachlässigt, jedoch nach ‘Ausführung der 
Integration in der ursprünglichen Form bestehen läßt, was ohne Beeinträchtigung der Ge- 
nauigkeit für die praktisch in Betracht kommenden Flugbahnhöhen möglich ist. Dadurch 
kann man nach Lösung der ersten Differentialgleichmg die Ordinate y durch eine einfache 
Quadratur als Funktion der Geschwindigkeit v erhalten: 


v 


® Wr ) 
= = — le . 0... ee AR 
, 15 + R er 
vo 
Ebenso gewinnt man dann die Zeit # durch eine einfache Quadratur: 
"d v 

u ı RU U 
» (18) 


vo 
Um schließlich noch die «-Ordinate zu gewinnen, geht man am zweckmäßigsten von 
der Komponentengleichung aus: 
dx dx 
di dv? 
de_v ı/[a\ (Wr N 
3: V (2) (+1) a a ey a 


woraus wiederum x durch eine Quadratur zu erhalten ist: 


=v-.cosd=v- Vı-( 





"rtP), 
g , 
oder 





v 


[2 .y(e)-rrı).ar. BEE 


vo 
Alle diese Differentialgleichungen gelten streng für ein beliebiges Widerstandsgesetz : 
Wrew,y)=ec: öly)v’- K(v,y). 


Will man jedoch die Auflösungen vorstehender Differentialgleichungen vollständig streng 
ohne die oben gekennzeichnete Vernachlässigung der zusätzlichen Höhenabhängigkeit von Wr 
und g durchführen, so muß man aus der zweiten Differentialgleichung (16) für y=dy:dv 
den Wert von p ausrechnen, was ohne Schwierigkeiten möglich ist, da diese Gleichung in p 
linear ist, den gewonnenen Wert von p nochmals nach ® differentiieren und die beiden so er- 
haltenen Ausdrücke für p und dp:dv in die erste Differentialgleichung (15) einführen, wo- 
durch diese in eine -Differentialgleichung übergeht, die ausschließlich die Variablen y und v 
enthält. Hat man daraus dann y als Funktion von v gewonnen, so kann man aus der zweiten 
Differentialgleichung für y durch einfache Substitution der gewonnenen Funktion y(v} auch 
die Funktion p(v) erhalten, so daß die Ermittlung der weiteren Ordinaten in der oben ge- 
kennzeichneten Form durch Quadraturen völlig allgemein und streng durchgeführt. werden 


kann. 
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$ 3. Vereinfachung und Umformung der Hauptgleichung. 


Da nun für eine Festlegung der Normalkomponente W, des Luftwiderstandes noch jede 
Grundlage fehlt, setzen wir zunächst in gebräuchlicher Weise in unserer neuen ballistischen 
Hauptgleichung W, =0. Wir setzen also damit voraus, daß die Symmetrieachse der Geschoß- 
form mit der Bahntangente stets zusammenfällt. Weiterhin kann man die neben der Einheit 
kleine Größe: 

9: R 


2 Saat vg 


stets durch einen konstanten Mittelwert ersetzen. Wir schreiben dann die Gl. (15) in folgender 
Form: 


„It W, r 
(g? > W = = +(2 . Wr 3.) .d® no 
A Br > 9 — ._ +3. ts. Wr+ wi. v—=0 (Ri). 








dv v p g 
Die Differentialgleichung hat somit die allgemeine Form: 
Ytha) y+hR@a) yet. 


Diese Gleichung würde ohne das von y freie Glied eine Bernoullische Differentialgleichung 
sein mit n=— 1, sie wäre also in y"’:2 linear. Wir formen daher unsere Differentialgleichung 
in folgender Weise durch Einführung neuer Variablen um, indem- wir setzen: 





dt de; Fr a? Eieg: dn=p-dp. 
Streng würde gelten: 
aut v_ von 
nn B-9 


N . i » a ) a „Ipe. 
Damit läßt sich schreiben, wenn wir obige Gleichung mit a dividieren: 


v? 


dn < Wer Wu 9, 72 m , Wr = v? 
dut?nt+ a; = — | y7-g+ (MT - 3.) ir. 
Mit: 
Se... 
—, dWr_dW + 
v=u.ct%,; em l4e) _— 
wird also: ; 
MWe_m 
1+8., ur 


somit endgültig: 


a 20 
dn 
nm] 
du g 1+ (23). 


En 


u 
(7 Wr-3,) 74° m 


Dadurch ist 7 als Funktion von « bestimmt und durch die Substitution: 





p=Vn; u=(l+2)-In— 
somit auch p als Funktion von u und von v. 


Wir können nun aber noch einen anderen schönen Weg zur Umformung unserer neuen 
Hauptgleichung einschlagen. Wir gehen dazu auf die obige allgemeine Form (14) derselben 
zurück und führen folgende Substitution ein: 


a EERREE (24). 
Man überzeugt sich dann leicht, daß in sehr einfacher Weise ganz streng und allgemein: 


pP + Wa=9-(g2+Wr)-2+ 9, v-1— 2) 
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wird. Somit erhält man folgende außerordentlich einfache Form der Hodographengleichung 
in den Variablen 2 und »: 


R _ 1 + Ey _ ® a Wx a 
-(9g2+ Wr) 2+ — + 9—Qgy:—| (1 — 2°) = . -y1-<2’; 
(25). 
BEER 
d+2)'97-90°7 I | 
Schließlich können wir noch statt der Variablen v» wiederum den Invo=uw als unab- 
hängige Variable einführen, wodurch sich die Form unserer Differentialgleichung noch weiter- 
hin vereinfacht: 





1 
—(gz+ Wr). 4, +9. 1-2)=Wy.yi-e rer 











Diese sehr einfache und übersichtliche Form der Hodographengleichung gilt also voll- 
kommen streng ohne jede Vernachlässigung auch unter Berücksichtigung der Höhenabhängigkeit 
von Luftwiderstand und Erdschwere nach einem beliebigen vorgegebenen Gesetz, trotzdem 
die partiellen Ableitungen des Widerstandes W7 nach % nicht explizit in der Gleichung auftreten. 

Als zweite Gleichung tritt hinzu eine Differentialgleichung für y: 


dy ® dy P 2 

dv pP’ du gz+tWr 
Wie man dieses System simultaner Differentialgleichungen in den Variablen 2, « und y ganz 
allgemein auflösen kann, wurde bereits oben im Anschluß an die Formen (15) und (16) der 
ballistischen Hauptgleichung dargelegt und ist ohne weiteres auf die vorstehenden Differential- 
gleichungen sinngemäß zu übertragen, wobei Wr und g vorgegebene beliebige Funktionen der 
Variablen y und « sind. 

Da die ausschließlich und allgemein in der Ballistik gebräuchlichen Methoden zur Er- 
mittlung des Luftwiderstandes Wr aus den erschossenen Flugbahnpunkten irgendeine.Normal- 
komponente W,„ vollkommen vernachlässigen, obschon diese auch bei drallstabilisierten Ge- 
schossen zumal bei größeren Abgangswinkeln erheblichere Werte erreichen kann, vor allem 
auch im Anfang der: Flugbahn, so müssen die auf diesem Wege erhaltenen Widerstands- 
funktionen Wr bzw. K(v) von den realen Gesetzmäßigkeiten für den tangentialen Widerstands- 
vektor mehr oder weniger abweichen und können somit keine zuverlässige Grundlage für die 
Aufstellung eines Naturgesetzes über den Luftwiderstand als Funktion der Geschwindigkeit 
abgeben. Nach meinen Erfahrungen zeigen auch die experimentellen Flugbahnbilder der Ge- 
schosse ganz besonders zu Beginn erheblich unregelmäßige Widerstandsfunktionen infolge der 
zunächst noch stärkeren Pendelung des Geschosses, wenn dieses das Geschützrohr verlassen 
hat. Dabei kann man sogar ganz einwandfrei feststellen, daß eine wesentliche Normalkompo- 
nente W, des Luftwiderstandes vorhanden war, da die ersten Flugbahnpunkte oft größere 
Höhenwerte y besitzen, als dem Abgangswinkel (#,) selbst ohne Tangentialkomponente des 
Widerstandes entsprechen können, indem nämlich ganz deutlich und einwandfrei für die 
ersten Flugbahnpunkte yZ>x-tg%, ist. 

Aus der zweiten unserer Ausgangsgleichung (1*) ersieht man sofort, daß unsere Sub- 
stitutionsvariable 2 nichts anderes wie der Sinus des Tangentenwinkels # ist: sind=z. 

Man kann daher die obige wichtige Form der ballistischen Hauptgleichung auch un- 
mittelbar aus den beiden Ausgangsgleichungen (1) und (2) für die tangentiale und normale 
Beschleunigungskomponente des Geschosses auf folgendem Wege gewinnen: 

Wir multiplizieren die zweite dieser Gl. (2) mit cos# und erhalten so: 

C 


v.cos®: re co®?d—+ Wy-cosd#, 


Br re 2) 


worin noch: 


- 
0) 


cos#-dd=d(sind)=dz; cos#=yl-—2’ 
ist, Aus der ersten Ausgangsgleichung (1) ergibt sich aber: 
dv dev 
di= - —; - =— ; 
9g-snd + Wr 9:2+ Wr 
Dieses in die vorstehende Gleichung eingesetzt unter Berücksichtigung von: dv/v=du ergibt 
sofort unsere obige einfache und doch völlig allgemeine ballistische Hauptgleichung (26). 
Setzen wir insbesondere wiederum Wy =0, lassen also die Symmetrieachse der Gesehoß- 
form auf dem ganzen Flugbahnwege mit der Bahntangente zusammenfallen, so nimmt unsere 
Hodographengleichung in völliger Allgemeinheit folgende Form an: 


(28). 
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ER, rn 
du 2+Wr:g’ 

worin: 





Wr=e: ö(y) fw,)=c:ö(ly)-v’-K(v, y), 
c=ballistischer Beiwert, ö=Luftdichte und K (v) = ballistische Beiwertsfunktion ist. 
Diese Gleichung ist mit der bekannten Form der Hodographengleichung, wie man sie 
in den ballistischen Lehrbüchern findet, identisch. 
Als zweite Bestimmungsgleichung tritt hinzu: 
Br 3 a ar Wr 


ee 





yo rl 





wobei W,: und g beliebig vorgegebene Funktionen der Variablen «=1n» und y sein können. 
Weiterhin kommen .dann noch folgende gleichfalls ganz allgemeine und sehr einfache Be- 
stimmungsgleichungen für die einzelnen Bahnelemente hinzu: 


dv e 
= —(g2e+W7r); di= _——- =v, de; 
EEE "TO Iiade IC Base 
(29). 


da=..yl-#-7; d=arcsin 2 


Damit sind sämtliche Bestimmungsstücke der Flugbahn streng dürch sehr einfache Aus- 
drücke dargestellt. 

Die obige Form der Hodographengleichung unter Berücksichtigung einer Normal- 
kompofente: 


(+7 . dz +8. yvi-=—(1-2)-du=0 


legt den Gedanken nahe, einen integrierenden Faktor zu suchen, der diesen Ausdruck in ein 
vollständiges Differential verwandelt, so daß dadurch die allgemeine Integration sofort durch- 
führbar wäre. Diese Aufgabe muß z. Zt. einer späteren ballistischen Forschung überlassen 
werden. 

Interessant ist jedoch noch folgende Bemerkung: Die obige Differentialgleichung (16) 
für y, die sofort auf eine Quadratur führt, stellt den Wert y als Differenz in der folgenden 
Form dar: . 


v=-\7- Wr 
=) g lpw 
vo 


Der erste Summand der rechten Seite dieser Gleichung ist dabei genau der Wert der Höhen- 
ordinate y für den luftleeren Raum: 
YWac.= (ve — v?):2g. 

Unsere Theorie ergibt also die y- Ordinate der Flugbahnpunkte in aller Strenge und Allgemein- 
heit durch eine Korrekturrechnung ’aus den Werten für den luftleeren Raum in einer gwissen 
Analogie und Fortsetzung zu dem bekannten Angelrutenprinzip. für die Konstruktion der Flug- 
bahnparabel im luftleeren Raum. 

Der Scheitelpunkt der Flugbahn wird insbesondere für 2=0, der Schnittpunkt mit der 
x-Achse für y=0, also für 





v 
vw—v? ve W: 
2 (30). 


ae ee 
tl ® r g pe) dv 





29 


erhalten. 


$ 4. Ermittlung der Widerstandskomponenten Wr und Wx. 


Unsere neue Hauptgleichung kann vor allem auch dazu dienen, um in einfachster Weise 
die Flugbahn zonenweise zu berechnen, indem man von Punkt zu Punkt der Kurve intervall- 
weise den Wert der Widerstandskomponenten Wr und Wy konstant setzt. In umgekehrter 
Reihenfolge stellt dasselbe Verfahren eine Methode zur einfachen und genauen: Bestimmung 
der Widerstandsfunktion aus den erschossenen Flugbahnpunkten dar. 

Wir vereinfachen zunächst die Schreibweise unserer beiden Hauptgleichungen noch da- 


durch, daß wir setzen: 
Wr _ g. MWn 
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wobei T und N dimensionslose Größen werden. Wir bezeichnen im folgenden diese Größen 7 
und N als „relative“ Widerstandskomponenten. Damit nehmen dann die Gl. (26) und (27) 
die Form an: 





n dz 2 7 / rear [9 
(+1). ,-U-)+N Vi-e=0 a re 
..: U. USE EP | > ORER: MOD. 
du g3s+T°' dw) &Bps+r (33). 











Die sehr einfache Form der zweiten der Gl. (33) weist darauf hin, daß es zweckmäßig sein 
wird, zur Bestimmung der Flugbalinordinaten als unabhängige Variable das Quadrat (v?) der 
Geschoßgeschwindigkeit zu wählen. Aus dieser Gleichung ergibt sich sofort der Wert T der 
relativen, tangentialen Widerstandskomponente: 





FEN 1 d (v?)\ 
1-43, 00). rd en re 











Sind nun die Flugbahnordinaten x, y und t durch photogrammetrische Flugbahnvermessung 

am geeignetsten mittels zerhackter Leuchtspur gegeben, so gelten Nie Beziehungen: 
; 4x a_Zer2ly 
I ee "= Ft 

Man bildet also zweckmäßig für das elementare Zeitintervall At der Vermessungspunkte 
aus den zugehörigen ersten Differenzen Ax und Ay nach vorstehenden Gleichungen die z- 
und v’-Werte und gleicht sie graphisch mit den Werten der y-Ordinaten aus. Der so durch 
eine Kurve festgelegte funktionale Zusammenhang zwischen 2 und »? einerseits und y und »? 
andererseits wird dann dazu benutzt, um nach Gl. (34) die relativen Werte T der tangentialen 
Widerstandskomponente nach bekannter Methode zu berechnen. . Die gewonnenen T- Werte 
kann man nochmals graphisch ausgleichen und erhält dadurch unmittelbar die relative Wider- 
standsfunktion T in Abhängigkeit von dem Quadrat v’ der Geschoßgeschwindigkeit. 


Nachdem so der T-Wert ermittelt ist, wird man die erste Hauptgleichung (32) dazu ver- 
wenden, um den relativen Wert N der hormalen Widerstandskomponente zu berechnen. Es 
ergibt sich unmittelbar aus Gl. (32): 


(35). 














no, a_ +8 u : 
N=y 1 z yi u du (36), 
wobei man die elementaren Differenzen Au aus der Beziehung: 
Au=lnv, ee ?x,+4y)-In(4’,;+4Y)) . -» -» -» . 89 


gewinnt. Es lassen sich also auch die N-Werte unmittelbar auf Grund des nach vorstehend 
beschriebener Methode experimentell ermittelten funktionalen Zusammenhanges zwischen 2 
und v? berechnen. 

Fällt die Symmetrieachse der Geschoßform mit der Bahntangente zusammen, so muß 
sich aus dieser Gleichung der Wert N=0 ergeben. Daher läßt sich dann T auch aus dieser 
Gl. (33) unmittelbar berechnen: 


du 


=). g5-=-(1+ 1 4) 


39 dy ur. 
Da der Normalkomponente im allgemeinen, solange N<1 ist, nur eine virtuelle und keine 
reelle Verschiebung in ihrer Richtung entspricht, so kann N ausschließlich eine Funktion des 
Anstellwinkels der Symmetrieachse des Geschosses sein. - Wenn also einmal die rein mathe- 
matische Theorie entwickelt sein wird, welche die Abhängigkeit zwischen N und a festlegt, 
so kann man damit nach vorstehender Methode aus den experimentell ermittelten N-Werten 
auch die Anstellwinkel a des Geschosses in jedem Flugbahnpunkt genau festlegen. 

Die Größe, die in der Aerodynamik schlechthin als Auftrieb bezeichnet wird, ist die in 
der Ördinatenrichtung senkrecht nach oben, also entgegengesetzt der Erdbeschleunigung, 
wirkende Kraftkomponente. Diese ist aus unserer Normalkomponente W, durch Multi- 
plikation mit cos® zu erhalten und also durch die Größe Wy :cos®# gegeben: Da cos # eine 
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Funktion der Geschwindigkeit » ist, wird somit der aerodynamische Auftrieb eine Funktion 
der Geschwindigkeit und des Anstellwinkels. Die vorstehende Beziehung lehrt ferner, daß 
dieser Auftrieb sich als das Produkt zweier Funktionen darstellen läßt, von denen die eine 
nur vom Anstellwinkel «a und die andere von der Geschwindigkeit (v) abhängt. 


$ 5. Ermittlung der Geschoßbahn durch eine neue mechanische 
Integriermaschine. 


Wir leiten noch einige wichtige und interessante Beziehungen für die Ordinaten der 
Flugbahn ab, indem wir das Geschwindigkeitsquadrat (v’) als unabhängige Variable einführen. 


Aus unserer Ausgangsgleichung (1): 
dv : } 
ZZ. Wr — 9: sindö=—g-(T-+2) 


gewinnen wir wie oben das Zeitdifferential: 


dv 


Dann lassen sich die beiden Geschwindigkeitskomponenten in den ÖOrdinatenrichtungen in 
folgender Weise umformen: 


up" ' 





dx — dx u 2 
gr wosd=v.yl—e; nr Far . (39); 
© > na u u 8 
gr sind=v:; dv gz+T Bi (40). 
Da aber: v-do=4di) ist, so folgt weiterhin: 
dx 1 ‚vi-®. dy _ 1 z dy_ 2 


(41). 








WG 24T’ dw) 2 + Ti|lde yi-e 





Die letzte dieser beiden Gleichungen ist mit Gl. (33) identisch. Wir setzen nun zur 
Vereinfachung der Schreibweise: 


29-.de=dx; 2g-dy=dy 
und können dann schreiben: 


u a (+) 


ea art Aw Hay 





(42). 


In analoger Weise formen wir die ballistische Hauptgleichung (32) um, indem: wir be- 
achten, daß: du=dv:v=d(v?):2»? ist: 





he ee, 
dw) z+1 z+-T 
Durch Einführung der Gl. (42) in den vorstehenden Ausdruck ergibt sich dann endgültig 
folgende wichtige Beziehung: 





de _„.dz 1-® dy 
de?) dw?) 2 dw) 








2v? A 











Die Gl. (42) und (43) können dazu’ dienen, um in einfacher Weise durch mechanische 
Integration mittels der von der Firma Askania in Berlin auf meine Veranlassung unter maß&- 


geblicher. Mitwirkung von R. Sauer’), Aachen gebauten Integriermaschine mit Hilfe der 


Variablen: v’, &, y und z=sin # die, Geschoßflugbalın unmittelbar aufzuzeichnen. Zu diesem 
Zweck formen wir die genannten Gleichungen in folgender Weise um: 


3) Vgl. R. Sauer und H. Pösch: Integriermaschine für gewöhnliche Differentialgleichungeu Z. VDI Bd. 87 
(1943), 8. 221. 
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Wr 
"=-2|| gdy+ | day) ESPRIT EEE . (44,); 
_[NM-» 
= \ 2 a ie RE er 


 Wı ig 1—# 
'L_ \ „ (dz—| .- Fe en 5 6 WE 
Darin wählt man zweckmäßig y als unabhängige Variable. Die ihm entsprechende 
Achse der Integriermaschine wird also unmittelbar und primär angetrieben. Man zeichnet 
dann x als Funktion von y mittels des Resultattriebes der Maschine auf. Desgleichen können 
mit weiteren Resultattrieben auch die Größen v* und z als Funktion von x unmittelbar auf- 
gezeichnet werden. Die empirischen Widerstandsfunktionen W7 (v?) und Wy: v?” werden durch 
Funktionstriebe in den Rechenvorgang eingeführt. Dabei ist aber Voraussetzung, daß die 
funktionalen Beziehungen: 
Wr yi-#. g1i—# 
ee 2 FE Zi 
durch Multiplikationstriebe gebildet werden 
Für die meisten praktischen Fälle wird man in den Gl. (44) die Erdbeschleunigung g 
wenigstens für den Integrationsvorgang konstant setzen können, so daß sich die Gleichungen 
durch folgende Schreibweise weiterhin vereinfachen lassen: 








Re - 
v”— — y+\-.49; Ba: ; 2 a ne 
= | HR .05; ZeB92.. 2222222. 4); 
j wi = ’ı11-e: ,_ 
= | Hd \ 3 ——-dy R . R e% : 5 rar (45,). 


. 


Unter der Voraussetzung von geeigneten Multiplikationstrieben benötigt man zur mecha- 
nischen Integration nach vorstehenden allgemeinen Gleichungen nur 4 Integrationstriebe, 
2 Funktionstriebe und 2 bzw. 4 Additionstriebe. Die funktionalen Beziehungen : 

2 2 2 
ı u -=-V()-1; Are. 
2 2 2 2 2 
kann man auch durch Kurven festlegen und mittels zusätzlicher Funktionstriebe in den 
Integrationsmechanismus einführen. Da zur Zeit brauchbare mechanische Multiplikatoren von 
der erforderlichen Genauigkeit nicht existieren, so kann man statt ihrer zusätzliche Inte- 
grationstriebe anwenden. 

Da die vorstehenden Gl. (45) noch relativ komplizierte Funktionen unter den Integral- 
zeichen enthalten, wollen wir noch einen anderen ebenso allgemeinen, dabei aber wesentlich 
einfacheren Weg zur mechanischen Integration der Flugbahngleichungen angeben. Wir ge- 
winnen diesen Weg leicht dadurch, daß wir eine zusätzliche Hilfsvariable: 


ne a 


t 
in die Gl. (45) einführen und der Größe z ihre ‚ursprüngliche Bedeutung als sind geben. Wir 
können die erforderlichen funktionalen Beziehungen aber auch unmittelbar aus unserer Aus- 
gangsgleichung (2) ableiten: 


dd 
v7 = Wyx — 09-008 v, 


worin wir unter Voraussetzung konstanter Erdbeschleunigung noch: 


. Se df 
2gdt=di; 2gi-f; dü-5 


schreiben können. Auf diese Weise erhalten wir: 


49-8. 
(2 


(N — ce0s9)-dt= (N — cosd)-du=(Wy —g-cosd)-du. . . . . (4). 
Hinzu nehmen wir noch die elementaren Beziehungen für das Bogenelement ds: 


dx dy dz dx dy (48) 
= s= DE ea ie \ 


—eos® sind 


ds 


” 7 e0s# sin®’ 
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wenn s=2yg:s ist. Endlich haben wir für die Geschwindigkeit (v) die Definitionsgleichung: 
— d 8 ” 
di r 


Dadurch haben wir folgendes System äußerst einfacher und doch strenger und all- 
gemeiner Gleichungen für den mechanischen Integrationsvorgang gewonnen. 


ds=v-dt=v’-du: ds=r-di=2vr-dü . .... 0.9). 


g = const g = variabel 
s=2-\v’dü ; = ) En ee a 9 Fe 
= \ (N — cosd) dü; 9= \(Wyr —g9:c08d):-du . .... (50,); 
®=—y—\Tds ; v”——2 \ gdy—2 \ Wr -.ds ... . &0,; 
y- | sind-ds s y= ee (50,); 
= | cosd.-ds I ea ar 


Als primäre Variable und unmittelbare Antriebachse der Maschine wählen wir die Größe u. 


Die Funktionen sind und cos®=sin (3 — 9) können wir durch eine Kurve mittels eines 


Funktionstriebes in den Integrationsmechanismus einführen. Da cos# aus sind durch eine 


Phasenverschiebung um — hervorgeht, so genügt an sich eine Sinus-Kurve, die durch 


. . . . . . . . ZT . 
einen Funktionstrieb mit doppelter Abtastung mit einer Verschiebung von d, => ein- 
geführt wird. Durch 3 Resultattriebe werden »v*, y und x als Funktionen von « aufgezeichnet. 
Die Kurve v? (u) wird sodann durch rechnerische Radizierung in die Kurve » (u) umgezeichnet. 
Diese Kurve wird in einem zweiten Integrationsvorgang mit nur einem Integrationstrieb dazu 
verwandt, um die Zeitvariable t als Funktion von « aufzuzeichnen gemäß der Beziehung: 


i=2-|v-da; t=|v-du ee a a rare A 

Da in Gl. (50,) für # auch in der rechten Kolonne die Erdbeschleunigung g in dem 
Ausdruck g-cosd neben Wy unbedenklich für fast alle praktischen Fälle konstant gesetzt, 
sowie in Gl. (50,) für v” das Integral \g:dy rechnerisch streng ausgewertet werden kann, 
wenn g nur Funktion von y ist, so benötigt man zur Durchführung der mechanischen Inte- 
gration insgesamt 5 Integrationstriebe, außerdem 2 Additionstriebe und 4 bzw. 6 Funktions- 


triebe, je 2 für die empirischen Funktionen T(v?) und N (v?) bzw. Wr(v’) und Wy (v*) und 
die theoretisch gegebenen Funktionen sind und cos# bzw. für die Funktionen g:cos® und 


2.|g-dy. 
Will man den Integrationsprozeß auch für variables g streng durchführen, so muß man 
die Gl. (50,) in folgender Weise zerlegen: 
= | Wy-du— | cosPdü, 
= | g-du. 


Man benötigt also 2 zusätzliche Integrationstriebe. Überdies kann man auch die Funktionen 
sind und cos# durch einen Integrationsmechanismus gewinnen, zu dem man nochmals zwei 
zusätzliche Integrationstriebe benötigen würde. Um schließlich die Höhenabhängigkeit der 
Widerstandsfunktion Wr entsprechend dem Ansatz: 


Wr (v’, y)= W, (y)- W, (ve?) 
zu berücksichtigen, ist noch eine Zerlegung des Integrals in Gl. (50,) erforderlich: 
| W,-ds=|\W,(r’)dw; w=|\W,(y)-ds. 


Dadurch würde sich dann die Gesamtzahl der Integrationstriebe auf 10 erhöhen, zu denen 
2 Additionstriebe sowie 4 Funktionstriebe für die Funktionen: 


we); WW; We’); gy) 


treten. Setzt man aber in der bisher üblichen Weise die normale Widerstandskomponente 
N=0, so reduzieren sich die erforderlichen Hilfstriebe für konstante Erdbeschleunigung auf 
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einen Additionstrieb und 3 Funktionstriebe. Aus den aufgezeichneten Kurven für %, & und £ 

gewinnt man die Raum-Zeitkoordinaten 4, x und t durch Division mit der doppelten Erd- 
beschleunigung entsprechend den Gleichungen: 

a nut u 

y=3,% —..s = 77: 


Dabei kann noch der jeweilige Wert von g entsprechend seiner Höhenabhängigkeit so- 





7 2 J 7 2 Jg + 17 
F-fv?da [rds ITfvl N=Nfvl sind cası8 Z«[cos Id3 
Se[N-cosH} du Iefsin das 

Bild 3. 


wie bei Fernbahnen auch gemäß seiner Änderung mit der geo- 

graphischen Breite gewählt werden. Damit sind dann sämtliche 

Variablen &, y, v sowie auch # und s als Funktion der Zeit t 

graphisch ermittelt. 

7 { .— Die Bilder 3 und 4 zeigen die zu dem mechanischen Inte- 
grationsvorgang gemäß den Gl. (50,) bis (50,) und (51) gehörigen 
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2 ER Schaltungsskizzen der Integriermaschine. Die Resultattriebe sind 
de v-v/u/ ER darin fortgelassen. Sie können beliebig an je 2 Achsen ange- 
Bild 4. schlossen werden. 536 


Eingegangen am 28. 3. 1944. 


Zur Theorie der Integriermaschine für gewöhnliche 
Differentialgleichungen. 


(Mitteilung aus dem Institut für Praktische Mathematik der‘ Technischen Hochschule Aachen.) 
Von R. Sauer und H. Pösch in Aachen. 


Für die auf dem Bush-Prineip beruhende neue Askania- Integrier- 
maschine wird die Frage nach dem Wirkungsbereich ‚durch Aufzeigung 
der sämtlichen mit dem Gerät lösbaren Differentialgleichungen vollständig 
beantwortet und eine Klassifikation der Schaltungen gegeben. 


8 1. Einleitung. 


Zur automatischen Lösung von Systemen gewöhnlicher Differentialgleichungen wurde in 
Anknüpfung an die amerikanische Bush - Maschine !) vor kurzem ein neues Universal-Integrier- 
gerät von den Askania-Werken in Berlin gebaut. Über die zahlreichen Anwendungs- 
möglichkeiten des Geräts in Mathematik und Technik und über die praktischen Fragen der 
Handhabung ist bereits an anderer Stelle”) berichtet worden und soll demnächst eine zu- 
sammenfassende Darstellung von H. Pösch erscheinen. Demgegenüber stehen in der vor- 
liegenden Untersuchung theoretische Gesichtspunkte im Vordergrund und zwar insbesondere 





1) V. Bush: Thedifferential analyzer — A new machine for solving differential equatione. J. Franklin Inst. Bd. 212 
(1931), S. 447/88. — S. Rosseland: Mechanische Integration von Differentialgleichungen. Naturwiss. Bd. 27 (1939), 
S. 729/35. — A. Walther: Neuzeitliche mathematische Maschinen. Elektrotechn. Z. Bd. 61 (1940), S.-33/36. — W. Meyer 
zur Capellen: Mathematische Instrumente. Leipzig 1941, S. 295/209. 


. 2) R. Sauer und H. Pösch: Integriermaschine für gewöhnliche Differentialgleichungen. Z. VDI Bd. 87 (1948), 
S. 221/224. 
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die Frage nach dem Wirkungsbereich der Integriermaschine. Sie wird durch Aufzeigung der 
sämtlichen mit dem Gerät lösbaren Differentialgleichungen vollständig beantwortet ($ 6). 
Vorher wird eine Prinziperläuterung des Gerätes ($ 2) und nach Einführung von Schaltplänen 
($ 3) und Kraftflußplänen ($ 4) eine Klassifikation der Schaltungen gegeben ($ 5). An ein- 
fachen Beispielen werden die Ausführungen verdeutlicht. 


$ 2. Prinziperläuterung der Integriermaschine. 


Die Integriermaschine ist nicht ein ein für allemal fest vorliegendes Gerät, sondern be- 
steht aus gewissen Grundtrieben (Integraltrieben, Summentrieben und Funktionstrieben), die 
in beliebiger Anzahl vorgegeben sind und bei jedem zu lösenden System von Differential- 
gleichungen in bestimmter Weise verkoppelt werden müssen („Baukastenprinzip“). 

Die unabhängige Veränderliche « und die abhängigen Veränderlichen %,, Y3:--., Yn 
werden dargestellt durch die Drehwinkel rotierender Achsen. Die «-Achse wird von außen 
durch einen Motor angetrieben und überträgt über die zwischengeschalteten Grundtriebe ihre 
Drehung zwangsläufig auf die übrigen Achsen derart, daß die Änderungen der Drehwinkel 
der sämtlichen miteinander verkoppelten Achsen gerade den vorgegebenen Differential- 
gleichungen entsprechen. 

Auf die konstruktive Verwirklichung’) im einzelnen, insbesondere auf die Ausstattung 
mit elektrischen Fernübertragungs- und Verstärkersystemen, brauchen wir hier nicht näher 
einzugehen. Es genügt folgende Prinziperläuterung der Grundtriebe (Bild 1): 














Funktionstrieb: v-f{u) 
„Jhtegraltrieb: rv- const [du 


Summentrieb: to-2m“us0 


Bild I. Prinzipskizze der Grundtriebe. 


1. Integraltriebe. 

Die Drehung einer Scheibe a überträgt sich auf ein scharfkantiges Reibrad b, so daß 
zwischen den Drehwinkein « der Scheibe und w des Reibrads und dem Abstand v des Reib- 
rads vom Scheibenmittelpunkt die Beziehung 


dw= const-v du 
besteht. Bei vorgegebener Drehung der Scheibe durch die «-Achse („Antrieb“) und vor- 


gegebener Verschiebung ® des Reibrades über eine Spindel von der v-Achse aus („Eingang“) 
ergibt sich für die w- Achse („Ausgang“) 


w = const - | vdu. 
Die Anfangsbedingungen der zu lösenden Differentialgleichungen werden an den Eingängen 
der Integraltriebe eingestellt. 


2. Summentriebe. 

Sie übertragen nach Art eines Ausgleichgetriebes eines Kraftwagens die Summe oder 
Differenz der Drehwinkel zweier Achsen auf eine dritte Achse, liefern also die Summe oder 
Differenz 

£ w=u+tv. 

3. Funktionstriebe. 

Die v- Achse dreht die Trommel ce, während die «- Achse über eine Spindel den Wagen d 
längs der Trommel verschiebt. Die eine der beiden Achsen wird aus dem Gerät heraus an- 
getrieben („Eingang“). Die andere Achse („Ausgang“) muß von Hand so nachgedreht werden, 
daß der mit Lupe ausgerüstete Wagen sich längs einer auf der Trommel vorgezeichneten 
Kurve bewegt. Die Drehwinkel der beiden Achsen sind dann durch die Beziehung 


v—f(u) 


verknüpft. Das Nachdrehen von Hand kann z. B. durch lichtelektrische Abtastung auto- 
matisiert werden. 
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Wenn man die Lupe durch eine Schreibvorrichtung ersetzt und beide Achsen aus dem 
Gerät heraus antreiben läßt, zeichnet der Wagen die Funktion v=f(«w) als Kurve auf. Auf 
diese Weise wird der Funktionstrieb zum Resultattrieb für die Entnahme der gesuchten 
Funktionen „;=f;(x). 

8 3. Schaltpläne. 


Bei Hinzunahme beliebiger Übersetzungen der Achsendrehungen und beliebiger Maßstab- 
faktoren der Drehwinkelmessung in den Veränderlichen x, y; erhält man nach $ 1 die Grund- 
beziehungen 


dw=kvdu, w=k|vdu EEE et 
5 ae ee ee a Me 
„wur: . :» in (3) 


für die drei Typen der Grundtriebe (Integraltriebe, Summentriebe, Funktionstriebe). 

Das zu lösende Differentialgleichungssystem muß in » Grundbeziehungen (1), (2), (8) 
zwischen (n-+1) Veränderlichen &, %,, Y2, --:, 4„n aufgespalten werden. Diesen » Grund- 
beziehungen sind dann ebenso’viele Grundtriebe des Geräts zuzuordnen, durch welche die den 
(n+1) Veränderlichen zugehörigen (n +1) Achsen zwangsläufig verkoppelt werden. 

Jedem lösbaren Differentialgleichungssystem entsprechen somit gewisse Schaltungen von 
(n-+1) Achsen über n Grundtriebe. Zur übersichtlichen Darstellung der Schaltungen kann 
man nach dem Vorgang von Bush Schaltpläne folgender Art verwenden, wobei für die Grun& 
triebe :die in Bild 2 angegebenen Symbole treten: 

1. Die (n-+1) Achsen werden als parallele waagrechte Linien übereinander gezeichnet 

mit Angabe der zugeordneten Veränderlichen x, %;. 

2. Darunter werden nebeneinander die Integraltriebe und Funktionstriebe gezeichnet, 

wobei die Eingangs-, Antriebs- und Ausgangsachsen als lotrechte Linien hinzukommen. 

3. Hierauf verknüpft man die waagrechten und lotrechten Achsenlinien entsprechend 

den vorliegenden, aus dem Differentialgleichungssystem gewonnenen Integral- und 
Funktionsbeziehungen. ä 
4. Zum Schluß werden die Summentriebe zur Verwirklichung der Summenbeziehungen 
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Bild 3 (rechts). Schaltpläne und 
Kraftflußpläne für „=f(x)-+g(y). 
d) Krafiflußpläne 








Diese Schaltpläne zeigen unmittelbar, wie die einzelnen Grundtriebe miteinander ver- 
bunden werden müssen, und dienen daher als Anweisungen zur Einstellung der Schaltungen. 
Sie lassen sich auch auf eine Schalttafel übertragen, indem durch Bedienen von Schaltern 
oder Steckern waagrechte und lotrechte Achsen der Schalttafel gekoppelt werden. 

Als Beispiel sind in Bild 3 unter (a) zwei Schaltpläne für die Differentialgleichung 

rn. 
angegeben. Beim linken Schaltplan ist (4) aufgespalten in 
(J):dp=fde, (J):dgq=gd«, 
(FR): f=f@), (F): g9=9W), 
(S):y=p+4, 
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beim rechten Schaltplan dagegen in 
(J):dy=hdx, (F):f=f@), (F)g=oWy), (S):h=f+9. 

Man sieht hieraus, daß einer Differentialgleichung mehrere Schaltungen entsprechen können, 
von denen bei der praktischen Anwendung des Gerätes die jeweils günstigste auszuwählen 
ist. Beim vorliegenden Beispiel ist der rechte Schaltplan der einfachere, da er einen Grund- 
trieb weniger als der linke enthält. Jedoch muß keineswegs stets die Schaltung mit der 
Mindestzahl von Grundtrieben die für die Praxis günstigste sein, sondern unter Umständen 
ist aus Maßstabs- oder technischen Gründen (z. B. Kraftverstärkung) eine andere Schaltung 
vorzuziehen. ö 


$ 4. Kraftflußpläne., 


Die Schaltpläne zeigen unmittelbar, in welcher Weise die Grundtriebe verkoppelt werden 
müssen. Ob jedoch die Schaltung richtig ist, d. h. ob. tatsächlich jede Achse zwangsläufig 
und eindeutig angetrieben wird, kann erst durch Verfolgung des Kraftflusses von der x-Achse 
über die Grundtriebe zu allen übrigen Achsen festgestellt werden. 

Der Kraftfluß läuft stets von den Antriebs- zu den Ausgangsachsen. Die Integraltriebe 
haben eine bestimmte Antriebs- und eine bestimmte Ausgangsachse. Die Antriebsachse treibt 
die Scheibe, die Ausgangsachse wird vom Reibrad getrieben, der Kraftfluß läuft also von 
der Scheibe zum Reibrad; eine Umkehrung des Kraftflusses könnte nur durch eine Anpressung 
des Reibrades erzwungen werden, die bei einem Präzisionsgerät naturgemäß unzulässig ist. 
Bei den Funktionstrieben kann der Kraftfluß sowohl von der Trommel zum Wagen laufen 
als auch umgekehrt, d. h. irgendeine der beiden Achsen wird aus dem Gerät heraus an- 
getrieben und ist Antriebsachse, die andere Achse wird von Hand oder lichtelektrisch be- 
tätigt und ist Ausgangsachse. Bei Verwendung eines Funktionstriebs als Resultataufzeichner 
werden beide Achsen aus dem Gerät heraus angetrieben und sind Antriebsachsen. Ähnlich 
sind bei den Summentrieben irgendwelche zwei der drei Achsen Antriebsachsen, die rest- 
liche dritte Achse ist Ausgangsachse. 

Die Eingangsachsen der Integraltriebe regeln die Verschiebung der Reibräder und be- 
einflussen hierbei den übrigen Kraftfluß nicht. 

Um den Kraftfluß unmittelbar ersichtlich zu machen, kann man die Schaltpläne in 
Kraftflußpläne transformieren, wobei es sich empfiehlt, für die Grundtriebe die in Bild 4 an- 
gegebenen Symbole zu benützen. Die den Kraftfluß führenden Achsen (x-Achse, Antriebs- 
und Ausgangsachsen) werden durch dicke Linien mit Pfeilen in Kraftflußrichtung dargestellt, 
die nur der Regelung dienenden Eingangsachsen der Integraltriebe durch dünn gestrichelte, 
ebenfalls bepfeilte Linien. Im einzelnen erfolgt die Konstruktion der Kraftflußpläne fol- 
gendermaßen: 

1. Zunächst wird die vom Motor angetriebene x-Achse als waagrechte Linie eingezeichnet. 

2. Hierauf zeichnet man die an die x&-Achse angeschlossenen Grundtriebe mit ihren 

Antriebs- und Ausgangsachsen ein und fügt Schritt für Schritt im Laufe des Kraft- 
flusses über die Antriebs- und Ausgangsachsen die weiteren Grundtriebe hinzu. Da- 
bei setzt man die hintereinander geschalteten Grundtriebe untereinander, die parallel- 
geschalteten nebeneinander. 

3. Zum Schluß werden die der Regelung dienenden, den übrigen Kraftfluß nicht stören- 

den Eingangsachsen der Integraltriebe durch dünn gestrichelte Linien dargestellt. 


Die Schaltung ist in Ordnung, wenn der Kraftfluß durch alle Grundtriebe genau einmal 
hindurchläuft und die Eingangsachsen aller Integraltriebe geregelt werden. In Bild 3 sind 
unter (b) die den Schaltplänen (a) entsprechenden Kraftflußpläne der Differentialgleichung (4) 
angegeben. 

Die hier zunächst für den praktischen Zweck der Verfolgung des Kraftflusses entwickelten 
Kraftflußpläne bilden ein sehr brauchbares Hilfsmittel auch für theoretische Untersuchungen. 
Mit ihrer Hilfe werden wir in den folgenden Abschnitten die Schaltungen klassifizieren und 
die Frage nach dem Wirkungsbereich der Integriermaschine, d. h. nach der Menge aller lös- 
baren Differentialgleichungssysteme, vollständig beantworten. 


$5. Klassifikation der Schaltungen. 


Die Schaltungen lassen sich klassifizieren nach der topologischen Struktur der Kraftfluß- 
pläne. Man erhält hierbei folgende Hauptklassen : 

1. Reihenschaltungen. 

Der Kraftfluß läuft in einem einzigen Zug durch alle Grundtriebe hindurch. Die Regelung 
der Eingänge der Integraltriebe erfolgt entweder aus dem vorangehenden Teil des Kraftflusses 
(Vorwärtskoppelung) oder aus dem nachfolgenden Teil (Rückkoppelung). 
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Bild 5 (rechts). Reihenschaltung, Parallel- 
schaltung, Verzweigungsschaltung. 


2. Parallelschaltungen. 

Der Kraftfluß gabelt sich an der «-Achse in mehrere parallele Zweige und durchläuft 
jeden dieser Zweige in Reihenschaltung. Die Eingänge der Integraltriebe werden entweder 
innerhalb der einzelnen Zweige geregelt (V.orwärts- bzw. Rückkoppelung) oder von einem 
Zweig zu einem anderen (Querkoppelung). 

3. Verzweigungsschaltungen. 

Im allgemeinen wird sich der Kraftfluß an mehreren Stellen gabeln und dadurch zu 
einer vielgestaltigen Verzweigung nach Art eines Stammbaumes führen. Hinsichtlich der 
Regelung der Eingänge der Integraltriebe kann man auch hier Vorwärts, Rück- und Quer- 
koppelungen unterscheiden. 

In Bild 5 ist je ein Beispiel einer Reihenschaltung, einer Parallelschaltung und einer 
Verzweigungsschaltung angegeben. 

Zu einer weiteren Unterteilung für jede der drei Hauptklassen kommt man durch die 
Unterscheidung in reduzible und irreduzible Schaltungen: 


a) Irreduzible Schaltungen. 

Eine Schaltung heißt irreduzibel, wenn kein einzelner Grundtrieb und kein Aggregat 
von Grundtrieben weggenommen werden kann, ohne daß der Kraftfluß oder die Regelung 
unterbrochen wird. Es ist also nicht möglich, einen Teil der Grundtriebe nach Abschaltung 
der anderen für sich allein laufen zu lassen. 


b) Reduzible Schaltungen. 

Eine Schaltung heißt reduzibel, wenn ein Grundtrieb oder ein Aggregat von Grund- 
trieben ohne Zerreißung des Kraftflusses oder der Regelung weggenommen werden kann. Eine 
reduzible Schaltung läßt sich demnach in Schaltungsteile folgender Art zerlegen: Ein oder 
mehrere Sehaltungsteile werden von der «-Achse angetrieben und sind im Kraftfluß und in 
der Regelung von den anderen Teilen unabhängig, können also jeder einzeln für sich in Be- 
trieb gesetzt werden. An diese Schaltungsteile sind die übrigen Schaltungsteile über einzelne 
Verbindungsachsen zur Kraftflußführung und zur Regelung. angeschlossen. Sie lassen sich 
ebenfalls einzeln für sich in Betrieb setzen, falls man die Verbindungsachsen durch geeignete 
Funktionstriebe koppelt. 





Bild 6 zeigt eine zweiteilige reduzible Schaltung. Der erste Schal- 
tungsteil ist an die x- Achse angeschlossen und entspricht dem Differential- 
pe" gleichungssystem 

U> dy, dp, _ 


3 2 % dz 9 dz 9ı- zer EUREN Fe ee TA 


Der zweite Schaltungsteil ist über die y,- und %,-Achse als Verbindungs- 
achsen an den ersten 'Teil angehängt und führt auf die Differential- 
Y gleichungen 





dı dı 
%—y, N—y,. er ee _ - 


Yu Uh 
dy, 


WE Yu TER dy, 
N“ ze. die sich nach Einsetzen der Lösungen y,=y,(x) und 9,=1,(x) von (5) 
chaltung. . . 
n integrieren lassen. 
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In reduziblen und in irreduziblen Schaltungen können Aggregate enthalten sein, die 
sich durch einen einzigen Funktionstrieb ersetzen lassen. In diesen Aggregaten kommt der 
Kraftfluß durch eine Verbindungsachse als Antriebsachse aus dem restlichen Schaltungsteil 
an und läuft in ihn durch eine weitere Verbindungsachse als Ausgangsachse zurück ; außer- 
dem erfolgt die Regelung ‚aller Integraltriebeingänge innerhalb des Aggregats. Wir bezeichnen 
derartige, mit Funktionstrieben äquivalente Aggregate als Funktionsaggregate. 

Naturgemäß entspricht jede irgendwie definierte Klasse von Schaltungen einer analogen 
Klasse von Differentialgleichungen. So sind beispielsweise den Schaltungen, in denen jeder 
Integraltrieb für sich allein oder zusammen mit einem- Funktionstrieb ein Funktionsaggregat 
bildet, durch Quadraturen lösbare Differentialgleichungen zugeordnet. 


86. Wirkungsbereich der Integriermaschine. 


Auf Grund der Kraftflußpläne kann man die Frage nach dem Wirkungsbereich der 
Integriermaschine durch Angabe aller lösbaren Differentialgleichungssysteme vollständig be- 
antworten. Man ordnet zu diesem Zweck jeder Ausgangsachse der Integraltriebe, Summen- 
triebe und Funktionstriebe eine Veränderliche &,, %,, ..., 4%. zu und stellt die aus dem Kraft- 
flußplan folgenden, durch die » Grundtriebe bestimmten Grundbeziehungen auf. 

Für die Integraltriebe liefert (1) unmittelbar Beziehungen der Form 


dy=k;ysdz oder dy=k;ady, oder dy=kydyı oder dy=krdz . (la). 
Bei den Summentrieben ergibt sich aus (2) durch Differentiation 
dy=adya+bidz ode dg=sadyatbdm : .» : 2.0.0. (2a) 
und ebenso aus (3) bei den Funktionstrieben 
dy=flaadz oder dy=flu)dyY - -» > > 2 222... (88) 


Indem man nun von einer Achse y; aus entgegengesetzt zum Kraftfluß bis zur «-Achse zu- 
rückgeht, erhält man eine Folge von Differentialbeziehungen (1a), (2a), (3a) zwischen dy;, 
dYk,...,d9Y,, dx. Nach Elimination von dyx,..., dy, bleibt eine Differentialgleichung von 
der Form 

dy=(YıtYot+t: + Ym)dx 


übrig, in der die Y;; Produkte von Funktionen (x), y(y,), -.-, 2(Yn) sind. Hieraus folgt 


Satz I: |} Jedes mit der Integriermaschine lösbare Problem läßt sich auf ein Differential- 
gleichungssystem von der Form 





Pn 


d pı d Da 4 
Yı _ > . Bu A Ü Yn _ > wi; en 
de , N Hat dz us, ee dz nd 3 “ z (‘) 


i=1 i=1 


& 














zurückführen, wobei die Summanden Y;; Produkte von Funktionen 9 (x), y, ,), 
..., X(Yn) sind. Die Zurückführung erfordert unter Umständen auch Diffe- 
rentiationen. 





Umgekehrt kann man zu jedem Differentialgleichungssystem (7) einen Schaltplan und 
Kraftflußplan angeben, d. h. in Umkehrung von Satz I gilt 


Satz II: | Jedes Differentialgleichungssystem (7) kann mit der Integriermaschine gelöst 


werden. 








Durch dieselben Schlüsse kommt man zu analogen Sätzen, wenn man nur gewisse von 
den drei Typen der Grundtriebe zuläßt. So folgt für Schaltungen ohne Summentriebe 


Satz III: [| Durch Schaltungen mit Integraltrieben und Funktionstrieben allein lassen sich 
die Differentialgleichungssysteme 





dy,__ 


dy, 
dz 


m Er dm x 
Y, dz 7 Dr ...y dx a. Y,. u a 3 ER, Zn (8) 











lösen, auf denen die rechten Seiten Y; Produkte von Funktionen @ (&), y (y,), - - -, 
x(y„) sind. Das System (7) spezialisiert sich zum System (8), wenn sämtliche 
Summen auf den rechten Seiten von (7) eingliedrig sind. 
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Für Schaltungen ohne Funktionstriebe erhält man den 


Satz IV: Durch Schaltungen mit Integraltrieben und Summentrieben allein lassen sich 
die Differentialgleichungssysteme 














dy, _» dp dm _p 

— . = u — ( 

2 2. 22 =P,, + = —E II SU 
lösen, bei denen die rechten Seiten P; Polynome in x, %,, -.-, %. sind. Läßt 


man auch keine Summentriebe sondern lediglich Integraltriebe zu, so sind alle 
Polynome P; eingliedrig. 





j » y dy; » 
Aus den n Gl. (9) und den aus einer der Gl. (9), etwa — P;, durch wiederholte 
KL 


Differentiation bis zur n-ten Ordnung sich ergebenden weiteren (n — 1) Beziehungen kann 
man alle Funktionen y, und alle Differentialquotienten PL für k=Fi eliminieren und 
erhält dann für y=y; eine Differentialgleichung »-ter Ordnung in der Form 

Pi@,4,9;:--, =; 


un 


hıerbei bedeutet P ein Polynom in x, 9, y', ..., Y Umgekehrt kann jede Differentialgleichung 
P=0 auf ein System (9) zurückgeführt werden. Zu diesem Zwecke setzt man etwa 


Gau FR a BR ee a a 
und erhält dann aus ; 
dP s 
Ge Pat Puy + Pyy’ ++ Pam yrt? = 0 
durch Auflösung nach y"”*” die zusätzliche Differentialgleichung 
TE Fr 53 
in der A und B wieder Polynome in &, %, %,,..., 4n bedeuten. Schließlich führt man mit 
Yn+.—=1/B 


noch eine weitere Veränderliche ein und spaltet (11) in die beiden Differentialgleichungen 


Ya Aynzı Ynsı= - By De u u Se u Dan, (12) 
auf. Hierin ist 

Bo Bx+ B,y, + Byı9: +. Byn Yn+ı 
ein Polynom in &, 9, %,, ---, Yn» Yn#+.. Somit ist das aus der Differentialgleichung P=0 her- 
geleitete System der (n-+2) Differentialgleichungen (10) und (12) mit den (na +2) gesuchten 
Funktionen %, %,, ---, Yn, Yn +. von der Form (9), w. z. bew. w. Wir fassen das Ergebnis zu- 
sammen in 


Satz V: || Eine Differentialgleichung n-ter Ordnung läßt sich dann ünd nur dann durch 
eine Schaltung der Integriermaschine mit Integraltrieben und Summentrieben 
allein lösen, wenn sig auf die Form 


|P@yy\,....y”)=0 er a 


gebracht werden kann, wobei P ein Polynom in x, 9, %', ..., 4°” bedeutet. 











Durch die vorangehenden Sätze ist die Menge der mit der Integriermaschine lösbaren 
Differentialgleichungen und Systeme von Differentialgleichungen vollständig abgegrenzt. Je- 
doch erfordert die Zurückführung einer vorgegebenen Differentialgleichung auf die kanonischen 
Formen (7), (8), (9) oder (13) unter Umständen mühsame Umformungsprozesse und gegebenen- 
falls mehrfache Differentiationen. Es ist daher keineswegs immer unmittelbar zu sehen, ob 
und wie die Zurückführung auf eine der kanonischen Forinen möglich ist. Tatsächlich ist der 
durch (7) und auch der durch (9) gekennzeichnete Wirkungsbereich so weitreichend, daß man 
im mathematischen und technischen Schrifttum kaum eine außerhalb des Wirkungsbereichs 
liegende Differentialgleichung findet. So lassen sich etwa alle Systeme von Differential- 
gleichungen, in denen die Differentialquotienten als rationale Ausdrücke von elementaren 
Funktionen, Besselfunktionen, Kugelfunktionen, hypergegmetrischen Funktionen usf. gegeben 
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sind, auf die kanonische Form (9) bringen. Der Wirkungsbereich wird nur dann über- 
schritten, wenn in den Differentialgleichungen empirische Funktionen von 2 und mehr Ver- 
änderlichen vorkommen. Aber auch in diesen Fällen kann man vielfach die Anwendung der 
Integriermaschine durch Zuhilfenahme von Netztafeln ermöglichen. 


Als einfaches Beispiel der Zurückführung auf die kanonische Form (9) behandeln wir 
die Differentialgleichung 


ren en + 
Die Funktion cosx wird als partikuläre Lösung y, von 


dy, dy, 7 
dx —N,, dx yes Yı 


angesetzt, so daß man als kanonisches System (9) die vier Gleichungen 


dy dy, dy. dy > 
h, ra "I Fr 7 de Urt 


de! 
erhält. Außerdam kann man (14) auch auf (13) zurückführen. Zu diesem Zweck eliminiert 
man aus (14) und der durch Ableitung gewonnenen Beziehung 


[77} 


y"=-sine+2yy 
die Veränderliche x und findet dabei 
Pa,y,Yy,y",y)=(y" —2yy)’ + WW" - PP’ —-1=0. 


Wenn die Zurückführung einer Differentialgleichung auf eine kanonische Form nur mit 
Hin?Zunahme von Differentiationen möglich ist, so liefern die zugehörigen Schaltungen die 
Lösungen der gegebenen Differentialgleichung nur bei Einstellung bestimmter Anfangs- 
bedingungen an den Eingangsachsen der Integraltriebe. Bei anderen Anfangsbedingungen 
ergeben sich die restlichen Lösungen einer übergeordneten Differentialgleichung, deren 
Lösungsmannigfaltigkeit die Lösungen der gegebenen Differentialgleichung als Teilmenge 
enthält. 


Der Wirkungsbereich der Integriermaschine wird nicht erweitert, wenn man als zusätz- 
liche Grundtriebe noch Produkttriebe zur Multiplikation von Achsen- 
drehungen hinzunimmt, da die Produktbildung auch mit Integraltrieben 
und einem Summentrieb durchgeführt werden kann. Um dies einzusehen, 
braucht man die Gleichung 

w=u-v 
nach Differentiation nur in die drei Beziehungen 


dp=udv, dgq=rdu, am=dp-+dyq 





aufzuspalten. Man erhält dann für die Produktbildung die in Bild 7 an- u RE 
gegebene Schaltung, bei der die «-Achse und die v»-Achse aus der übrigen „4; en 
Schaltung heraus angetrieben sind. schaltung. 


$ 7. Zusammenfassung. 


Die zunächst der Untersuchung des Kraftflusses in den Schaltungen der Integrier- 
maschine dienenden Kraftflußpläne führen auch zu theoretischen Erkenntnissen. Mit ihrer 
Hilfe wurden Klasseneinteilungen der Schaltungen vorgenommen und kanonische Formen ent- 
wickelt, auf die sich alle mit der Integriermaschine lösbaren Differentialgleichungssysteme 
zurückführen lassen. Dabei ergibt sich ein so weitreichender Wirkungsbereich, daß die meisten 
in der Mathematik untersuchten und in den Anwendungen auftretenden Differentialgleichungen 
in ihm enthalten sind. 519 


Eingegangen am.3. 11. 1943. 
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Eine Bemerkung zur numerischen Differentiation durch 
Approximation, ausgeführt am Beispiel der Kugelfunktionen 
als Approximationsfunktionen. 


Von Erwin Fehlberg in Berlin-Charlottenburg. 


Zur Berechnung der Koeffizienten, die bei der Entwicklung gegebener 
Funktionen nach Kugelfunktionen auftreten, wird die zu entwickelnde 
Funktion abschnittsweise durch Parabeln ersetzt und das Produkt dieser 
Näherungsparabeln mit den entsprechenden Kugelfunktionen exakt inte- 
griert. Das Verfahren bewährt sich auch bei der Berechnung eines 
Näherungswertes der ersten Ableitung der gegebenen Funktion, wie an 
einem Beispiel gezeigt wird. 


I. Allgemeines. 


Bei der Approximation einer stetigen Funktion f(r) durch Kugelfunktionen P, (r) erhält 
man die Koeffizienten der Entwicklung bekanntlich durch die folgenden Integrale 


1 
2 1 
FOR...» \ fd Pu()dr, BERN N. er rt 
1 


5) ’ 


Die Genauigkeit der Approximation hängt außer von der Gliederzahl im wesentlichen von 
der Genauigkeit ab, mit der die obigen Koeffizienten berechnet werden können. In der ein- 
schlägigen Literatur wird zur Berechnung dieser Integrale meist die Simpsonsche Regel 
oder dergleichen empfohlen. 

Während für kleine Indizes n die Genauigkeit dieser Verfahren bei praktisch uurch- 
führbarer Schrittweite der Integration im allgemeinen ausreichend ist, trifft dies bei größeren 
Indizes » nicht mehr zu. Der Grund hierfür liegt in der mit zunehmendem Index n 
sich ständig vergrößernden Schwingungszahl der Kugelfunktionen für das Grundintervall 
—1<rs-+l. 

Um bei der Berechnung der Integrale (1) die Schwingungen der Approximationsfunktionen 
ausreichend genau zu erfassen, wurde hier ein neuer Weg 
eingeschlagen: Das Grundintervall -1<r=s-+1 wird in 
Teilintervalle der Intervallbreite Ir=0,1 eingeteilt. Für 
je 2 aufeinander folgende Teilintervalle wird die Funktion 
f(r) durch eine Parabel 2. Ordnung ersetzt, die in den 
3 Endpunkten der Intervalle mit den Funktionswerten f (r,.) 
übereinstimmt, Bild 1. Für solche Parabeln sind aber die 
Integrale (1) exakt ausführbar. | | 






1 


| 
| 
| 
| 
1) 
} 





Der wesentliche Unterschied gegenüber allen früheren 3 
Verfahren besteht also darin, daß nur für f(r) ein an- 7° Tr Teer Une 
genäherter, für die Kugelfunktionen dagegen der genaue Bild 1. 


Funktionsverlauf benutzt wird. 

Durch Differentiation des Kugelfunktionsansatzes erhält man für den gesuchten Diffe- 
rentialquotienten von f(r) einen Näherungswert. 

An einem Beispiel wird gezeigt, daß die hier vorgeschlagene Metlıode zur Berechnung 
der Integrale (1) vollkommen ausreichende Näherungswerte für die erste Ableitung von f(r) 
liefert, während die Auswertung nach der Simpsonschen Regel (bei gleicher Schrittweite) 
für die Differentiation zu bedeutend schlechteren Ergebnissen führt. 

Der hier für die Berechnung der Kugelfunktionskoeffizienten eingeschlagene Weg läßt 
sich in gleicher Weise ganz allgemein für die Berechnung der Koeffizienten einer Entwicklung 
nach Funktionen eines vorgelegten Orthogonalsystems anwenden '). 


II. Die Berechnung der Kugelfunktionsk effizienten. 


Nach der Newtonschen Interpolationsformel erhält man für die Parabel 2. Ordnung, 
die durch die drei Punkte 


1) Dies wird zur Zeit von Herrn E.-A. Behrens für die Entwicklung periodischer Funktionen in Fourier- 
sche Reihen durchgeführt. Die Ergebnisse werden demnächst veröffentlicht. 
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T% ’ fir) ‚ft 
%ı, =. +dr ’ HJ =fkı=fet+ Afk; 
%Y=rt2dt, fe )ehanmht2Aft+ 4 fe 


ö A e 
ho=f+ en +zL ke — u) —%,,)- 


Einsetzen in (1) liefert: 


Tk+2 "k+2 Tk+2 


geht: 


2 Ei, [ 4 [ A 
nr ri > If \ Puär+ SR. \ (177) Pul)dr+ I. Me \ (T—T4)(t- tu) Pulnldrl. 
%k 





k Tk 173 

Setzt man abkürzend: 

Tk+2 

\PA.Odı=J, (rn, tı)=(- nr | r, (d)drz, 

ck T_(k+2) 

"k+2 "—k 

| @—u)Pu(d)dr= Kr, (u, = ("+ | ed )Pu(lo)dr, 

Tk T_(k+2N 

Th+2 T_k 

"a—rr)(e — 7%, ) Pal) dr= Ip, kr un (IR, tar) =D" | du) (T—T_ rn) Pule)dr 

°k T_(k +2) 


so findet man wegen 
Ps Ye [Pas.(0 - Ba (2)] . 


1 1 
In lin, m )=5, 41 [Pas (tr) — Pur (il 5, »>1 -[Pu-, (Tk+2) Pa-,(t%)] 


5) 
Krın On te), ri [Pu+ı (Tk4 2) — Pa (Tk+2)] 


1 
"Bn+l Uns: (Un, ter) 2; In-ı (tk, Tk42)] 
2(Ar) 
Ir, k+ı,n (Um Tkı 9) [Pas re) — Pn-ı (Ik+2)] 
4(Art) 1 1 
a rang [Pu sltePaltne] gu [Palin Pas (ixn)] 
Ar 


+5 21 nr (on, TRr) — In-ı (Tr, Trr>)] 





2 1 
+5, +1 mu PAR OR Bas A (Tr, Tk43)] 
1 


”- z Un (TR, tk) — J n— RR) 





(4), 


(5), 


Die Zahlenwerte dieser Funktionen sind den Zahlentafeln 1 bis 3?) zu entnehmen. Die 


Werte dieser Tabellen lassen sich aus jeder Tabelle der Kugelfunktionen’°) ableiten. 
Für den .Kugelfunktionskoeffizienten a, ergibt sich dann 
<0,1 


,=5 > {ir + (—- 2" f-,] an +1) In (tk, Trr:) 
k 


1 
rm [4 kr +(- 1" 4 len +1) Kr.n (ik, Tas.) 





1 
MPIPETL [If + —V" A’f-)2n +1) Lr,k+ı,n (Tks 1) 


2) Die Berechnung der Tabellen sowie die Zahlenreehnung zu Kap. III verdanke ich Frl. M. Kristahn. 
s) Z.B.: PEreNE: Tafeln der Kugelfunktionen P„(x) und ihrer abgeleiteten Funktionen. Acta Soc. Se. 
Fennicae 32 (1906), N. 6. Die Tabellen bei Tallquist gehen bis n=8; für die weiteren Werte bis n—12 wurden 
Tabellen für P„(x), Pan (x) und Pf (x) auf 5 Dezimalen mit einem Argumentschritt von 0,1 berechnet. Abzüge werden 
vom Verfasser gern zur Verfügung gestellt. 


(8). 
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’ Zahlentafel 1. 
2n+ 1) In (irrt) = (— D"(2n +1) Ju (t_ik+ 2), TR). 
(TR; tk+2) 
(0; 0,2) (0,2; 0,4) (0,4: 0.6) (0,6; 0,8) (0,8; 1) 
n 
0 + 0,200000000 + 0,200000000 - 0,200000000 + 0,200000000 + 0,200000000 
1 + 0,060000000 + 0,180000000 + 0,300000000  0,420000000 + 0,540000000 
2 — 0,480000000 0,360000000 0,120000000 + 0,240000000 0,720000000 
3 — 0,203000000 — 0,525000000 0.595000000 0,245000000 + 0,693000000 
4 + 0,587520000 + 0,123120000 0,503280000 0,686880000 + 0,479520000 
5 + 0,374924000 + 0,718212000 + 0.174460000 0,738892000 + 0,158796000 
6 0,601036800 + 0.315806400 + 0,760468800 - 0,315369600 0.159868800 
7 0,544926300 — 0,600646500 + 0,599878500 + 0,334897500 — 0,375140700 
8 + 0,539473920 — 0,713647080 — 0,194528280 + 0,796207920 0.427506480 
9 + 0,688168076 + 0,195201538 0.819840310 r 0,773187042 0,317185096 
10 — 0,42027246 + 0,85458908 0.62540633 + 0,29144786 0,10035815 
11 — 0,78588862 + 0,31685807 + 0,19159946 0,32833678 + 0,13408818 
12 + 0,2640153 — 0,6597033 + 0,7948116 0.6973162 + 0,2981926 
Zahlentafel 2. 
(Zn +DKrın (TR, )=(2Rn +1) (— Dr! Kr,n(T-(krds TR) 
(tk; tk +2) 
(0; 0,2) (0,2; 0,4) (0,4: 0,6) (0,6; 0,8) (0,8; 1) 
n 
0 + 0,020000000 | + 0,020000000 -+ 0,020000000 + 0,020000000 + 0.020000000 
1 + 0,008000000 + 0,020000000 + 0,032000000 + 0,044000000 + 0,056000000 
2 — 0,047000000 - 0,033000000 0,007000000 — 0,031000000 + 0,081000000 
3 0,026880000 — 0,056280000 0,057680000 0,014280000 + 0,090720000 
4 + 0,054420000 + 0,001836000 0.059388000 0,063732000 + 0,084564000 
5 + 0,048998400 + 0,070514400 + 0,001601600 0,083820000 + 0,065577600 
6 0,049794940 + 0,048289020 + 0,070915780 0.057524740 + 0,039445380 
7 - 0,069765120 0,045655020 + 0,077989080 + 0,003933180 + 0,012849840 
8 + 0,035346992 0.085151497 + .0,009401102 + 0,067133741 0.008315428 
9 + 0,085474099 0,009906285 0,072861975 + 0,096394310 0.020368480 
10 - 0,013666049 + 0,083360652 0,089697842 + 0,075041198 0,022523873 
11 0,093415372 + 0.067612681 0,022723998 + 0,014953480 0.016706663 
12 0,01198478 0,03945981 + 0,06628978 - 0,05098563 — 0,00662900 
Zahlentafel 3, 
(Zn+1) Ir,ksın (iR, TR) (2n +1) (— N” Ip,krın (T-(Kr Ds TR). 
(Tk;Tk+2) 
(0; 0,2) (0.2; 0,4) (0,4; 0,6) (0,6; 0,8) (0,8; 1) 
n 
0 | + 0,000666667 + 0,000666667 + 0,000666667 + 0,000666667 -+ 0,000666667 
1 | + 0,000400000 + 0,000800000 + 0,001200000 + 0.001600000 + 0,002000000 
2 .- 0,001486667 — 0,000886667 + 0,000113333 + 0,001513333 + 0,003313333 
3 — 0,001325333 - 0,002100000 — 0,001754667 + 0,000270667 + 0,004536000 
4 + 0,001470000 -— 0,000658800 0.002538000 0.001647600 + 0,005540400 
5 + 0,002351213 + 0,002152920 0.001063773 0,003282547 + 0.006165720 
6 0,000851162 + 0,002670018 + 0,001810614 0,003617198 + 0,006265350 
7 - 0,003202270 — 0,000373007 + 0,003614928 0.002135737 + 0,005761822 
8 — 0.000255032 — 0,003496346 + 0,002380058 + 0,000764540 + 0,004688905 
9 + 0.003658061 — 0,002412770 0,.001305217 + 0.003799279 + 0.003203693 
10 + 0,001676056 + 0,002114470 0,004367486 + 0.005356674 + 0,0015610% 
HM - 0,003569566 + 0,004440592 - 0,003803490 r 0,004366756 + 0.000055856 
12 0,003187451 -+ 0,001079996 + 0,000378518 + 0,001007242 — 0,001051262 
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Hierbei bedeutet f_, den Funktionswert f(,, an der negativen Stelle 
z_r und Af-x=f-w+u—f-r sowie I’f-»=Af-iwrv —Af-x- 


Unter Benutzung der Tabellen 1, 2 und 3 kann die Berechnung der Kugelfunktionskoeffizienten 
bequem mit der Rechenmaschine durchgeführt werden. 


III. Fehlerbetrachtung. 


Man hat zwei verschiedene Fehler zu berücksichtigen: 

1. Da man die vorgelegte Funktion f(r) auf dem Intervall <—1,+1) für die Berechnung 
der Kugelfunktionskoeffizienten a, durch ein System von aneinander gehängten Parabeln 
2. Ordnung f, (r) ersetzt, so ist: 

+1 


EEE 5 5 te 
—1 


2. Wenn man die Kugelfunktionsentwicklung der aneinandergehängten Parabeln mit 
dem Index n = N abbricht, so ist: 


+1 N ; 

ER h)—2 Et tere A 
—1 ed 
+4 N 
: u 

) — 2 er = | 

2, Im dr = gern, A, VRR ARE TE SE ar Sr ie (11). 
—1 n—=(0 


Die Konvergenz der Kugelfunktionsentwicklung — nicht nur im Mittel, sondern an jeder 
Stelle x des Intervalls <—1,-+1> — ist bekanntlich immer gewährleistet, wenn die zu ent- 
wickelnde Funktion der Dirichletschen Bedingung genügt, d. h. wenn sich das Intervall 
<—1,-+1> in endiich viele Teilintervalle so zerlegen läßt, daß die Funktion auf jedem Teil- 
intervall monoton ist*). Diese Bedingung erfüllt der Parabelzug f,(r) sicher, womit die Kon- 
vergenz der Kugelfunktionsentwicklung sichergestellt ist. 


Die Ausdrücke für 2, und 2, sollen noch weiter ausgeführt werden. Um einen 
Näherungsausdruck für 2, zu gewinnen, genügt es offenbar für f(r) eine Newtonsche 
Interpolationsparabel 3. Ordnung anzusetzen. Und zwar wähle man auf <r;, 7;,,> diejenige, 
die an den 4 Stellen r;._,, 7x, Tgy,, 7x4, mit f(r) übereinstimmt, während auf <r;,,,Txy4s> 
diejenige anzusetzen ist, die an den 4 Stellen 7x, Txy,, Tk4+3, 7%4, mit f(r) übereinstimmt: 


fr 1° fr 


A fi A fa. 
<Tp, dt ri ge De Pre en I Tkzı)(T-Te;,), 


If 1" fı f 


P ER , ir k 
(Fk kr: fl) Fat At er tar, (7; rer rg „erlernte, de Ter+2)- 


Einsetzen in (9) liefert: 


. "+ 
er A? ; 2 [ . r > 
k 2 
a Tk+2 
ENTE 
(81 -(Ar \ 1 


Tky+ı 


| 


- u’ (Tr, — TR dr, i 


Eine elementare Zwischenrechnung liefert schließlich: 


8 1 ” 
= . -At- fi Ve N PU 2). 
9-5 gi 4t 2 Bf... . (12) 
q 
Der Summationsbuchstabe i soll andeuten, daß jetzt über jeden Teilpunkt des Intervalls 
<—1,+ 1> zu summieren ist (bei dem Summationsbuchstaben k war nur über jeden 2. Teil- 


punkt zu summieren). 


4) Vgl. z. B. Frank-v.Mises: Die Differential- und Integralgleichungen der Mechanik und Physik, Bd. |, 
Braunschweig 1930, S. 438. 
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Für 2, erhält man, wenn man für f,(r) den Ausdruck (2) einsetzt nach kurzer Rechnung: 


92-4] N Inr+2- I Inl-un+g DInten+z LS unlten 


k k k k 
N ; (13). 
4 5” 2 ı 2 \ k 2 a! \ a - - 2 
+3 / Unf+p 2, enM- TE 
k k "n 0 


Aus dem Wert für 2, kann man nun offenbar sofort eine Aussage darüber machen, 
bis zu welchem Index N man die Kugelfunktionsentwicklung höchstens zu erstrecken braucht. 


Da nämlich 2, der Fehler ist, der dadurch entsteht, daß man die Kurve f(r) durch an- 
einandergehängte Parabeln 2. Ordnung f,(r) ersetzt, so hat es offenbar keinen Zweck, die 
Kugelfunktionsentwicklung von f(r) oder genauer des Parabelzuges f,(r) weiter fortzusetzen, 
als bis 2, etwa gleich 2, wird, weil man sich doch letzten Endes für die Entwicklung von 
f(r) und nicht für die des Parabelzuges f, (r) interessiert. 


Eine analoge Fehlerbetrachtung läßt sich nun für die Differentialquotienten von f(r) 
anstellen. Es seien für den ersten Differentialquotienten die Ergebnisse kurz mitgeteilt. 


La 
In 


ET 
ı 
Für @, ergibt sich ein etwas komplizierterer Ausdruck, da die ersten Differential- 
quotienten P/(r) der Kugelfunktionen nicht mehr orthogonal sind: 





= eh + gr Sant +g, Damen 
"5 0 Ui ET 
k k h 
N 
2 "s 2 3 ” Er2 Be 
u. ° \ a I a (4) It 4 fi) ® 2 P.trr)] 1 fı 48 fj|\ 
Ir er Like u un l -_ J 
== I k k 
N (15). 
2 wu " = \ 
. m ® 2 - ® 2 eo ws) +9 
+8 2 ns (If) 2n +1): Iunlta, ı% >)| 
n—=] k 
N’ 
+2 £ Pi Am n Am Ay + 2 A, n" a), 
1=sm,n=N n—| 
m>n 
+ 1 
Für die Zahlen Ayun= | Pm(r): P,(r)dr ergibt sich: 
BE 
I(l!+1), wenn m und n beide gerade oder beide on 
f sind und / die kleinere der beiden Zahlen ist. 16) 
ee | 0 , wenn von m und r die eine Zahl gerade, die andere | En  . 
Zahl ungerade ist. 
N 
Man kann die Entwicklung a, P,(r) offenbar wieder abbrechen, sobald 2; etwa 
nel 
gleich 2/ wird. Es ist natürlich nicht notwendig N’=N, d.h. man muß im allgemeinen die 
Kugelfunktionsentwicklung a, P„(r) für f(r) und die Reihe Na, P,(r) für f’(r) bei ver- 
n—=1 n=]1 


schiedenen Indizes abbrechen. 


IV. Ein Beispiel. 


Um den Unterschied zwischen den nach der Simpsonschen Regel und den über die 
Hilfsfunktionen J, K, I, berechneten Kugelfunktionskoeffizienten zu zeigen, genügt es, als 
Beispiel für f(r) eine Kugelfunktion selbst zu nehmen. Es wurde f(=P,(r) gewählt. 


Die Entwicklung ergibt, wenn bei beiden Integrationsverfahren die gleiche Intervall- 
breite Ar=0,1 zugrunde gelegt wird: 
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nach Simpson über J, K, L 


nach Simpson! über J, K, L 





Ay Ay Ay Ay 
0 + 0,00069 + 0,00069 7 N) 0 
1 0 0 8 + 0,39330 + 0,01309 
2 —+ 0,00889 +0, 0789 9 0 0 
3 0 0 10 + 0,70867 — 0,01251 
4 + 0,04648 +0,02219 11 0 0 
5 0 0 12 +0,91486 — 0,04044 
6 + 1,15957 + 1,02402 





Hält man daneben die exakten Werte a„=0 für n=+-6 und a,=1, so erkennt man 
deutlich, wie mit zunehmendem Index n die Ergebnisse der Simpsonschen Regel immer 


unzulänglicher werden. 


Die Fehlerbetrachtung für f()=P,(r) und für f()=P;(r) ergibt: 








2, = 0,0074 2! = 0,78 
sowie N 2, und N’ 2, 
6 ' 0,00039 6 0,82 
8 | 0,00037 8 0,93 
10 | 0,00036 10 | 0,76 
12 | 0,00023 12 0,10 


Daraus folgt, daß man die Entwicklung für P,(r) bereits mit dem Index N=6 abzubreehen 
hat; für P,(r) kann man entweder bis N'=6 oder bis N’=10 gehen. 











ar - 
00 er nach Sımpson berechnet IN’=&) 
-— = über Kl berechner(W’-6) 
—— nach Simpson berechner(N-6) Fe 
006-  ——- überJ,Kl berechnet(N:6) u 
u = ee I 70 
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Bild 2. Bild 3. 


In Bild 2 ist die Differenz AP, zwischen den genauen Werten und den nach beiden 
Methoden gewonnenen Näherungswerten P, dargestellt. 


In Bild 3 ist die entsprechende Differenz A P; aufgetragen. 


Der Vergleich der Differenzen der Abbildungen läßt die Reichweite des hier entwickelten 
Verfahrens erkennen. 533 


Eingegangen am 25. 1.44. 
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Bemerkung zu einem Satze von Hurwitz. 
Von Herbert Bilharz in Braunschweig. 
Als Erweiterung eines Satzes von A. Hurwitz wird im Anschluß an I. Schur 
ein Determinantenkriterium angegeben, welches notwendig und hinreichend 
dafür ist, daß sämtliche Wurzeln einer algebraischen Gleichung mit kom- 


plexen Koeffizienten nur negative Realteile besitzen. Zur numerischen Be- 
rechnung der einzelnen Determinanten wird ein einfaches Schema mitgeteilt. 


Man verdankt A. Hurwitz') den nachstehenden 
Satz 1: Die sämtlichen Wurzeln der reellen algebraischen Gleichung n-ten Grades 
n 
fi (&) 2a,0=0. TEN er 


‚—0 


mit a,>0 haben dann und nur dann negative Realteile, wenn alle » Hauptunterdeterminanten 


B: een 
der quadratischen Matrix 
a, a, DS 8 
a, A, 0 0 00 
— a a a, 0, a, a, 0 
M = 5 4 3 2 1 0 (3) 
Asn-ı Asn- ; E £ 5 ; ; Ay/ 


positiv sind; in M ist dabei a.=0 zu setzen für u>n. 


Unter Zuhilfenahme einer von E. J. Routh?) angegebenen Methode, der sogenannten 
kreuzweisen Multiplikation, lassen sich die Determinanten (2) mit geringem Aufwand nume- 
risch berechnen. Man geht hierfür von den beiden ersten Spalten der Matrix (3) aus und 
erhält durch sukzessive Bildung der nach höchstens n +1 Schritten abbrechenden Folge von 
einzeiligen Matrizen 


a 
-7 le, 0, 0, ...)=[0, b,, d,, du, -. .]; 


a 
[@,, &, @, -- J-75 [1 u... Je, 
1 





(4), 
b 
[B,, ba, du...) [Cı, 6 &, ...]=[0,d,d,,d,,...]; 
1 
[e; ; Ey, O5, +» 1-14, d,,d,, z J=1[9, e,, @3, 85, ic . 
die D, als Produkte: 
D,=a, D,=a,b,, D,=a,b,c,, D,=a,b,c,d,, er N ce 


Beweis: R bezeichne die Substitutionsmatrix ?) 


es ee ci 
R ne Bat a 
iz wre" 


1) A.Hurwitz: Über die Bedingungen, unter welchen eine Gleichung nur Wurzein mit negativen reellen Teilen 
besitzt. Math. Ann. Bd. 46 (1895), S. 273 bis 284. 

2) E. J. Routh: Die Dynamik der Systeme starrer Körper, Bd. 2. Leipzig (1898), S. 222 bis 234. 

3) Im Hinblick auf das Folgende ist der obige, an eine frühere Note des Verfassers (Geometrische Darstellung 
eines Satzes von Hurwitz für Frequenzgleichungen fünften und sechsten Grades, Z. angew. Math. Mech. Bd. 21 (1941), 
S. 96 bis 102) sich anschließende Beweis gewählt worden. Vermittels elementarer Determinanteneigenschaften ließe 
sich auch direkt zeigen, daß die wiederholte Anwendung des Algorithmus (4) di» Determinanten (2) in Dreiecksgestalt 
überführt. 
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mit 

b,=8,;, für v=0 mod. 2 

0<rv<n 
0,0, +, — Gyr 2 5, i u; 2 
= tr 7??? für»=1 mod.2 (7); 
a, | 
b, =0 


die Ausübung dieser Substitution auf D, liefert 


.. 7 
R(D,) = —- ee tes 5 MR 
D, 


Für v<2 ist (d,,) wegen 


D=a, D,=DR(OD)=aRla)=a,b, \, (9) 
erfüllt. Durch vollständige Induktion bleibt also zu zeigen, daß 
D,,,=a,b,c,...9,h, ee 
wird, wenn für »=1,2,..,n—1 
ee aan gar: SR 
das Produkt der durch (4) gegebenen » Faktoren a,, b,, €,, -.., 9, ist. Weil 
Rla)=b, Rb)=RR a) = b1 Be ;. Baba et, u 
Alc)=d,,..., R(g,)= h, | 
und 2 TIL DERN. L r 
Ria,b,e.....- = Rla) Rl5)-Rke)-... -Ag)=b,c.d... gb. . . . (1) 
gilt, folgt aus (8) die Behauptung: 
Bear 3 


Seit Erscheinen der Hurwitzschen Abhandlung sind für Satz 1 noch weitere Beweise 
veröffentlicht worden; unter ihnen verdient die originelle und im Gegensatz zu Hurwitz 
rein arithmetisch-algebraische Beweisführung von I. Schur‘) besonderes Interesse, da in 
ihren Hilfssätzen die Koeffizienten von (1) dem Körper der komplexen Zahlen angehören 
können und die sonst üblichen Hilfsmittel wie Cauchysche Indextheorie, quadratische 
Formen und Sturmsches Theorem nicht benötigt werden. 

In der vorliegenden kleinen Mitteilung soll unter der Voraussetzung komplexer 
Koeffizienten in (1) — wie sie u. a. bei Stabilitätsuntersuchungen linearer Schwingungs- 
systeme mit phasenverschobenen Kräften auftreten — ein Satz 1 entsprechendes Determinanten- 
kriterium angegeben werden, auf welches der oben beschriebene. Algorithmus zur numerischen 
Berechnung anwendbar bleibt. 

Zur Abkürzung heiße eine algebraische Gleichung (mit reellen oder komplexen Koeffi- 
zienten), deren sämtliche Wurzeln negative Realteile besitzen, eine Hurwitzsche Gleichung, 
oder auch kurz eine H-Gleichung. 

Dem Polynom n-ten Grades 


KhKaazu, ta, +a,2°+--. +2"... 2.2.20... 0, .,(18) 
mit komplexen Koeffizienten a, =0, 1,..., n) und a,=+0°) werde das Polynom 
=, - +4,20” 3. +H- Ya... 202200. (16) 


zugeordnet, wobei &,=a,—ia, die zu a,=a,-ia,, konjugiert komplexe Größe bezeichne ; 
ferner sei £,„ eine beliebige komplexe Zahl mit Re&,<0 und 


HH) act, ya) ta, dt :. : . (MM. 
Dann gilt nach Schur®) 
Satz 2: Es ist n 
N 2 W" 2°T Fe 


=d 


*) I. Sehur: Über algebraische Gleichungen, die nur Wurzeln mit negativen Realteilen besitzen. Z. angew. 
Math. Mech. Bd. 1 (1921), S. 307 bis 311. 

5) Wir sehen von einer Normierung des Polynoms (15) ab, um die Homogenität in den spater auftretenden 
Formeln nicht zu verletzen. 

6) Vgl. S. 310 in loe. eit. ®). 
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dann und nur dann eine H-Gleichung, wenn 


+0, Re (4) >0 TE a u 3 
a, 
und die Gleichung (r — 1)-ten Grades 
fo-i (X) =’ [fa On fi Yn] b, 7 b,a+b,x” +: +b, En d., Fasn » . (20) 


eine H-Gleichung ist. 
Die Bedingungen (19) lassen sich wegen 


a 1 1 
Re (2 =; Rea,ä,= ke 
a, la,/? | +a,4,) 
in die Ungleichung 
m BL FÜ N 


zusammenfassen. 
Mit Satz 2 ist ein rekursives Verfahren gegeben, das nach höchstens (n — 1)-maliger 
Anwendung zum Abschluß kommt. Wird nämlich Satz 2 auf Gl. (20) angewandt, worin 


nunmehr RE ö b 
rn -,)=b,—b,2+b,2°” 7... +1" b._,r te (3) 


und bei beliebigem &,_, mit Re&„_,<0O 
mn =bebn sr ctbin-n Yen mbacrbin-,atbin-, - - (24) 
gesetzt werden, so ist (20) genau dann eine H-Gleichung, wenn 
a ee in 


und 
fin (8) x" [Ifn-ı Pa-ı — fn-ı Yn-ı] =o+6, z<+,2%°+.+ Cn „an? ==0 u (26) 


eine H-Gleichung ist. Auf diese Weise erhält man für (18) als Verallgemeinerung der be- 
kannten: Stabilitätskriterien von Routh eine Folge von n Ungleichungen: 


a +ra4,>0, 55,+55,>0, Sa +anc>0, ...:.2 2.2.2.2, 


welche notwendig und hinreichend sind, damit die vorgelegte Gl. (18) eine H-Gleichung sei. 
Für die Koeffizienten in (20) findet man 


I: ” 2 Mä 7 m - Yı 
b,o+1: Wat — Mt + En tl r+ ti) — (A,e+ a +, ; 


n — 1]\ 
FF NE EIER = (28) 
\ | 2 | 
B ii ie _ 2 ' 
b30 = (Ad, 41 + A, 4 Ft Entldofsc+: ü A,d,0+3) (4, Q,04 1 Ad,Gyo +1)3 
und daraus für die entsprechenden Koeffizienten des Polynoms (23): 
r Be 2 ER : 2 | P 
b,0+1 = (0427 Q,d,0+3) — End, ,045 + A Agar 5) — (A, Ayo54+2 7 per >39 
> = 29). 
m (dä 7 £ Slä 7 A 7 
b,o = (00,04, + %lyar 1) — Enid, A,o+2 - A,Älyg+2) (a, Aya4ı £ a,@,041)} 
Die Ungleichung (25) lautet hiernach 
bb, +b,b, = — (Ent En) [(G,a, —a,d,)”’+(d,a, + a,ä,)(ä,a,+a,d,)— (d,a,+a,d,)}] >0 (30), 


wobei über £, mit Re&„<0 noch beliebig verfügt werden darf. Wählen wir &, negativ reell, 
ERBETEN, . . 3 0 rt e 
und setzen bei Berücksichtigung von (22) in (28) und (29) 
-1 


na Be ne ee a 
2(ä,a,+a,ä,) 


Fr 


so folgt aus (30) 
(d,a, — a,ä4,)’+(@,0,+0,4,)$(d,a,+a,4,) — (d,a,+a,G,)) 
(ä,a, + a,ä,) 


Analog erhält man aus den Koeffizienten von (26) für die Bedingung (27,) 
+ El, + End, b,—b, 5, +(b,b,+b,6,) fd, b,+b,5,)— (b,b, + b,b,)Y] >0 (84), 


worin mit Berücksichtigung von (33) 


bb, +b,b, >0. .. (8). 
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Er ’ 
gewählt werden kann, so daß 
_ - (b, b, Fe b, b,)? + (b,b, + b.b,)X(b, b, + b, b,) Er (b, b, + b, b,)} (Qp 
= a a aA ey Ta Ta 7 02/7 0. . (86). 
Mn (2b, + 6.6) Pr 
Die sukzessive Bildung der weiteren Ungleichungen in (27) ist nach den in (33) und (36) 
angeschriebenen Ausdrücken evident. Bezeichnet nämlich R diejenige Substitution, welche 
den Koeffizienten a, und ä, die Größen b, und b, gemäß (28), (29) und (32) zuordnet, so 
folgen die iinken Seiten in (27) durch wiederholte Ausübung von R auf (22): 


Rasa; +a,ä)=b,b,+bb, RRAa+aA)=Rbb +6) = +, ... (BT). 


An Stelle von (22) kann offenbar die Determinantenbedingung 


| a —iü,| 
D=| . Pe EA er 
—ia, ä, 
treten; wegen (38) und 
a, —iG, 0 0 | 
P.2%5 ” BIN Bi = en 2 
D=| va, for ER (0, 40," ai r E u | (39) 
| —% a, va, a,| +(4,a,+0,4,){(@,a,+a,&,) — (d,a,+ a,ä,)} 
| ii, —-—ü —a, iä, 
wird 
_ | 
65,+6d=T :- rennen. (ON) 
und daher (33) äquivalent zu 
Fe 534 


Um zu einem allgemeinen Bildungsgesetz zu gelangen, betrachten wir die 2n-reihige 
quadratische Matrix 


a,.— ii, 0 0 0 0 0 0 - 
@ ia, ä, a, —iä, 0 N) 0 0 - . \ 
— 4, vä, —ia, a, a, —iü, 0 0 
ia u a IR —-is ä a —IG, 
a, —tä, sa —ü, -—a, ia, —ia, ä, 
ee 0 2 (#2), 
-0 ti, —ia, ä, a, ia, ia, —äü, 
ia, —äü, —a, iä, —ia, ä, a, —iä, 
a, —iäü, ia —ü, —a, sa, —ta, dä, 








\ r ; - ; . ; ; i us 


worin a,—0 zu setzen ist für» >n. Von unwesentlichen Faktoren abgesehen, wird ihre Deter- 


minante D, gleich der Resultanten der beiden Polynome f(ix) und —if(-ix)=-if*(ie). 
Die Folge der Hauptunterdeterminanten aus (42) sei mit 
ae 


bezeichnet; für reelle a, sind diese Determinanten mit den Hurwitzschen D, aus (2) durch 
die Identitäten 


D,=%a,D,_,D,, D=(-2ifa.Dd (=l..,n-l) .... (4) 
verknüpft. 
Nach (40) wird 
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und analog (8) gili allgemein für v„=1,2,...,n—1 die elementar, aber nicht kurz zu be- 
weisende Relation 
D 
R(D)=7 (46). 


1 


Die linken Seiten der Ungleichungen (27) werden, falls noch D,=1 gesetzt wird, die 
Quotienten zweier Determinanten: 


D, 


% (47), 


D = ı _D mn . 
d, +04, =5 b, b,+bb,=- z 66, +%C,= 


0 ü D, . 
und wir erhalten in Verbindung mit (46) durch vollständige Induktion als Verallgemeinerung 
des Hurwitzschen Satzes den 


Satz 3: Die algebraische Gleichung n-ten Grades (18) mit beliebigen komplexen 
Koeffizienten a, und a,=F0 ist dann und nur dann eine H-Gleichung, wenn alle Haupt- 
unterdeterminanten gerader Ordnung 


in re ee MR 
der Matrix M positiv sind, 

Dieser Satz liefert wie Satz 1 n Determinantenbedingungen’). Im Vergleich zu den 
Hurwitzschen D, hat sich dabei die Ordnung der D, entsprechend den zwei Kompo- 
nenten a,, a, der Koeffizienten a, verdoppelt. Die Determinanten D, in der Folge (43) sind 
im allgemeinen komplexwertig und von Null verschieden; da Satz 3 über sie keine Aussage 
enthält, kann ihre Berechnung unterbleiben. 

Das für die Berechnung der D, aus Satz 1 gültige Rechenschema ist, wie wir nun- 
mehr zeigen wollen, auf unsere D, übertragbar. Durch einfache Umformung erhält nämlich 


D, die Gestalt 


a. -a 0 0 0 0 
ıi:E & u —a 0 0 
I— a, a’ a, a‘ u —a 
— u. —a —a, a, a” a‘ 
ı da -d -d -a —a a, 
Dd,=9.| & 4 « -a -a -a 2’ .d, (49), 
— 46, a, a, a, a, —a, 
a —ad —a, a, a, a, 





XUM 


wonach die Hauptunterdeterminanten gerader Ordnung der Matrix M sich von den entsprechen- 
den Hauptunterdeterminanten gerader Ordnung der Determinante d, nur durch positive Fak- 
toren unterscheiden, welche bei anderer Wahl der £, aus Satz 2 nicht aufgetreten wären und 
im folgenden beiseite gelassen werden können. Beiläufig sei bemerkt, daß zwischen den D, 
und den Hauptunterdeterminanten ungerader Ordnung aus d„, welche; wir in Analogie zu. (43) 
mit d,, d,, ..., d„_, bezeichnen wollen, keine derartige Beziehung besteht. 

Um festzustellen, ob eine vorgelegte Gleichung den Bedingungen aus Satz 3 genügt, ist 
ohne Beschränkung der Allgemeinheit also notwendig und hinreichend, daß alle n Deter- 


minanten gerader Ordnung 
er en a ai 


positiv werden. Weil diese Determinanten prinzipiell den gleichen Aufbau wie die D, aus (2) 
haben, wird die eingangs mitgeteilte Rechenvorschrift anwendbar, welche sich folgendermaßen 
aussprechen läßt: 

Die von oben nach unten gelesenen Elemente der zweiten Spalte in d, werden von links 
nach rechts als Horizontalreihe geschrieben und darunter als zweite Zeile die entsprechenden 
Elemente der ersten Spalte gesetzt. Die Elemente dieser zweiten Zeile multipliziere man mit 


”) Mit Hilfe des Sturmschen Satzes ließe sich überdies zeigen, daß die Anzahl der Zeichenwechsel in der Folge 
D,(>0), Da, ..., Du für eine beliebige komplexe Gleichung n-ten Grades gleich der Anzahl der Wurzeln mit positivem 
Realteil ist. Nach eventueller Elementartransformation der unabhängigen Veränderlichen könnten die D, also auch 
zur Bestimmung der Anzahl der in einer von Rex<0 verschiedenen Halbebene oder der in einem Quadranten der kom- 
plexen Ebene gelegenen Wurzeln verwandt werden. 'Ein diesbezügliches determinantenfreies Kriterium für die 
Halbebene /mx>>0 wurde kürzlich von H. Wielandt angegeben. 

6 
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einer Zahl m, derart, daß naeh gliedweiser Addition dieser Produkte zu den darüberstehenden 
Zahlen deP ersten Zeile die Summe der ersten Glieder Null wird. Diese wird fortgelassen 
und die Summen aus den anderen Elementen, um einen Platz nach links versetzt, als dritte 
Zeile angeschrieben. Mit der zweiten und dritten Zeile dieses Schemas wird nun das Ver- 
fahren ‚wiederholt: man multipliziere die Werte der dritten Zeile mit einer Zahl m,, so daß 
bei gliedweiser Addition dieser Produkte zu den darüberstehenden Zahlen der zweiten Zeile 
die Summe der ersten Glieder Null wird. Das Ergebnis aus den anderen Elementen wird — 
ohne die Null an erster Stelle — um einen Platz nach links versetzt, als vierte Zeile auf- 
geschrieben, usw. Spätestens mit der (2n»+2)-ten Zeile ist das Verfahren beendet; ein vor- 
zeitiges Abbrechen kommt bei H-Gleichungen nicht vor, und bei ausnahmsweisem Auftreten 
in anderen Gleichungen läßt es sich durch geringe Abänderung der numerisch gegebenen 
"Koeffizienten vermeiden?). 

Die aus diesem Algorithmus resultierenden Elemente der ersten Spalte seien mit ö,, Ö,, 
u, 8, ), EN Ön-ı, On, u =0 bezeichnet. Die Hauptunterdeterminanten von d, folgen 
daraus durch fortgesetzte Produktbildung: 


d, => Ö,, d, = d, Ö, o d, == d, Ö, r d, zur d, Ö, z d,  — d, Ö, ” | 
. x ‚ j . (51). 
A En; 4,=4 Ö,, o., d,, BE RER RR d,= dy_, On 


Für eine H-Gleichung der Form (18) müssen nach (50) alle » unterstrichenen Ausdrücke 
positiv werden. 
Abschließend folge hierzu ein Beispiel: Gegeben sei die Gleichung fünften Grades 
f)=r’+ (10+ 392. +(4 +19)’ + (114 +50) 
+A1+Bd)r +10 40) 0 
Das oben beschriebene Rechenverfahren liefert für diese vorgelegte Gleichung das Schema: 
“- 20 171 45 - 44 —83 1 


(52). 











My, 

140 73 — 114 — 19 10 0 0,1428 
181,429 28714 — 46,714 — 1571 1 — 0,7717 
50,843 -77,953 — 17,788 9,228 — 3,5684 
306,883 16,761 — 34,500 1 — 0,1657 
— 80,730 — 12,072 9,062 3,8014 
— 29,199 — 0,052 1 — 2.7715 
— 11,928 6,291 — 2,4421 
— 15,415 1 —. 0,7738 
5,517 2,7941 
1 — 5,5170 

0 

Mit den Werten der ersten Spalte 
| ö, ö, Ö, d, Ö, Ö, Ö, ö, Ö, 6, 








— 20| 140 | 181,429 | 50,843 | 306,883 | — 80,730 — 29,129 — 11,928 — 15,415 5,517, 1 | 0 





erhält man nach (51) 
d,>0, d,>0, d,>0, d,>0, d,<0, d,>0, d,<®, d,>0, d,>0, d,>0, 
wonach (52) als H-Gleichung erkannt ist; ihre Wurzeln z,=x,+ix, sind 
BI LE IT FT 1 8 





1-1 —ı1| 1) — ii —5 


1 1 2| —2 4 0 
haben also in der Tat nur negative Realteile. 525 
8) Für den Anlauf der numerischen Rechnung ist zu beachten, daß in Gl. (18) — nach evtl. Multiplikation mit 


einer von Null verschiedenen (komplexen) Konstanten — die Voraussetzung ao == 0 durch ao’ =F U ersetzt werden kann. 


Eingegangen am 7. 12. 1943. 
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KLEINE MITTEILUNGEN 


Allgemeine Zerlegung der Beschleuni- 
gung des komplan bewegten starren 
ebenen Systems. 1. Aufgabenstellung. 
Es sind verschiedene spezielle Zertegungen des Be- 
schleunigungsvektors 5b , eines beliebigen System 
punktes A des allgemein bewegten ebenen starren 
Systems angegeben worden, von ge die von 
A.Schell!) und M.Grübler?) am bekanntesten 
sind. Die erstere wird durch die Gleichung d =, 

+ b,, die letztere durch 4 ,=b,-—+ b., *) dargestellt, 
worin mit # als Drehpol und r als "Tangentialpol 
b, AP) w - AT)e, während mit W als 
Wendepol % AP)e und bu=(AW)o? @& be- 
deutet, wie üb li ch, die Winkelgeschwindigkeit, 
e=o die Winkelbeschleunigung des komplan be- 
wegten ebenen Systems. Ö, ist unter dem Namen 
„Triebbeschleunigung”, b,, als „Wendebeschleuni- 
ung“ bekannt. Beiden Zerlegungsartei ist offenbar 
cemeinsam, daß die eine der beiden Beschleunigungs 
komponenten ais relative Normalbeschleunigung 
des Systempunktes 4 gegen einen anderen System- 
punkt (P bei 5, bzw. W bei db), die andere als 
relative Tangentialbeschleunigung gegen einen 
weiteren Systempunkt (7 bei db, bzw. P bei b, 
gedeutet werden kann. Diese Zerlegungsart ist 
einer weitzehenden Verallgemeinerung fähig derart. 
daß sich zu einem beliebigen Systempunkt M 
stets ein weiterer Systempunkt N so angeben 
läßt, daß bi nm ttan mit nun = (AM) w*, 
tın=(4AN)e nach der Bezeichnungsweise von 
F. Wittenbauer?°) ist. Es soll der Zweck dieser 
Arbeit sein, die Zuordnung der Punkte N zu den 
Punkten M anzugeben. 

2. Synthetische Begründung. In Bild 1 
bedeutet neben den bereits erwähnten Polen P.W, 
T, noch @ den Beschleunigungspol und A einen 
beliebigen Systempunkt, dessen Beschleunigung 
b4—= 44, mit AQ den „Beschleunigungswinkel“ 
einschließt, wobei {y=<{[ PTW ist. Zur Auf- 
tragung aller Längen sei in Pr 1 ein Längen- 
maßstab (1em der Zeichnung = a” der Natur) und 
zur Auftragung aller ee ein Be- 
schleunigungsmaßstab ( 1em der Zeichnung = yu/sek® 
der Natur) angenommen.*) Dann ist 








— {N} — € 
AQı =(AQ)— und AR,=(AQ 
Ist nun ein beliebiger Systempunkt M gewählt und 
durch Q, die zu QM parallele Gerade Q,M, bis zum 
Sehnitte M, mit AM gezogen, so wird durch AM, 
im gewählten Beschleunigungsmaßstabe y offe nbar 








1) W. Schell: Uber den Beschleunigungszustand des 
ebenen unveränderlichen, in der Ebene beweglichen Sy- 
stems, Schlömilchs Z. Math. Phys. Bd. 19 41871), S. 185; 
ferner W. Schell: Theorie der Bewegung und der Kräfte, 
1. Bd. 2. Aufl. Leipzig (1879), S. 441 bis 474 

2, M. Grübler: en l. Aufl. 
Ss. 1236 ff 

3) R.Beyer: Tec nalen he Kinematik, Leipzig 1931, S.240, 
241. Man vgl. auch _ ittenbauer: Graphische Dy- 
namik, Berlin 1923, 8 8, 39. Die Bezeie hnungsweise ist 
zunächst der des Bahn von R. Beyer angeschlossen. 

4, Es werden stets die natürlichen Längen in Klammern, 
z. B. (AP) mM, ihre zeichnerischen Bilder aber durch Quer 
striche, z. B. AP zur Darstellung gebracht, so daß also 

APym — AP.am ist. AP bedeutet also die dimensionslose 
Anzahl der em, in der in der Zeichnung mittels des Längen- 
maßstabes «aM die wirkliche Länge (AP)M zur Darstellung 


Berlin 1917, 


gelangt. Analoges gilt für die Beschleunigungen: Es ist z.B. 
b"; sek? — (AQ)M Vor+ „2sek £ AQ am Vos+e2 „sek? 


Daher ist das graphische Bild 
1Aı = ‚= AQ . aM 


» 


4=: „m/sek? Y .,a 7 ‚2sek — 
somit, wie es sein muß, eine dimensionslose Zahl. Hin 
sichtlich der Maßstäbe vgl. man Abschnitt ?. 





qN 
Lv 
RX , 
9% 
/. / / 
/ 1 
us / 
| / 
| / 
ab. 
/ 
/ 
pi 1} 
1 
72 
x RER - 
yM 


Bild 1. Zum Beweis der Zerlegungsmöglichkeit des Be- 
sehleunigungsvektorsb in» 4M +t 4N und Konstruktion 
der Punkte N zu vorge + benen Punkten M. 


die Beschleunigung (AM) o? dargestellt. die als die 
eine der eingangs erwähnten Besehleunigungs- 
komponent nten aufgefaßt werden soll. Da 5b,= 
AA, =AM,- -M,A,, so ist M,4, nach Richtung, 
Sinn und Größe im gewählten Beschleunigungs- 
maßstab >) offenbar die zweite Komponente, die 
im Systempunkt A aufgetragen, durch die Strecke 
AN, dargestellt wird. Soll diese Komponente nun 
als relative Tangentialbeschleunigung gegen einen 
gg Systempunkt N u, werden, so 
muß N. auf einem Halbstrahl AN _ı AN, offenbar 
in = her Entfernung von A an werden 
daß von ihm aus gesehen die Strecke AN, unter 
dem sogenannten „Winkelbeschleunigungswinkel“ # 
und zwar im selben Sinne drehend erscheint. wie 


IR,=(4Q) 


um den Beschleunigungspol @ dreht. In Bild 1 
fällt dieser Drehsinn mit dem Uhrzeigersinn zu- 
sammen, weshan .. N auf dem nach oben gerichteten 
Halbstrahle AN liegen muß. Es ist nämlich tg # 
der Winkelbeschleunigung e direkt proportional 
und der Proportionalitätsfaktor ist hierbei nur 
von den Maßstabskonstanten a” und ymisek? ab- 








die Beschleunigungskomponente 


hängig. Aus dem AQAR, folgt nämlich sofort: 
4 ( m 

ted = Au, _ Ad)e:r __« . sek? (1). 
AQ (AQ):a ym/sck? 


Ferner ist bekanntlich die Beschleunigung bp des 
Drehpohles P, die im Wendedurchmesser d liegt 
und die Richtung PW besitzt, durch die Beziehung 
bp=do?=ee gegeben, woraus 

[DJ e e 


= 22 =—=p . 6) 
r q 3 © DE Fu ı 





folgt, worin wieder die überstrichenen Größen sich 
auf die Strecken der Zeichnung beziehen. 
6*+ 
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Ferner ist AA.Q,M, »4AQN, denn es ist 





0, 4,=AR,=AQtgP# [Gl.(1)] und ıAQ 


(in Bild 1 zweimal gestrichen) und 


M,4,=4AN,=AN tgd und LAN 


(in Bild 1 einmal gestrichen), woraus zunächst folgt: 


<IQ, AM.=<[QAN 


Ferner muß auch M.Q, = N@tgö und LNQ sein 
(in Bild 1 gekreuzt), also auch MQ_LNG. Ferner 
ist, da 





und AM, =(AM) — : 


auch 


2m @* u w® 
M,Q,=(MQ)—- =MQu-. 
y y 
Somit folgt durch Gleichsetzung der beiden für 
M,Q, gewonnenen Ausdrücke und wegen (1) 





MQ, = Ma”- = Mte9—NQ - 8. 


/ 


Daraus folgt nun bei Beachtung von (2): 


——— u PR FE 
NQ=M9—-=MO- = Mlcoty, | 
€ d (3). 


NQ _ıLMQ im Gegensinne von & 


Man vgl. Bild 1, wo<{y=MNRQ bei N zur Kon- 
trolle eingezeichnet worden ist. Man kann die 
durch (3) angegebene Zuordnung der Punkte N zu 
den gewählten Punkten M so formulieren: 


Die Punkte N werden aus den Punkten 
M durch eine einfache Drehstreckung 
um den Beschleunigungspol@ gewonnen, 
wobeis die Drehung um 90° im Gegen- 
sinne der Winkelbeschleunigung e, die 
Streckung aber im Verhältnis des Glei- 
chenkreisdurchmessers e zum Wende- 
kreisdurchmesser d erfolgt. 


3. Vektorieller Nachweis. Dieser gestaltet 
sich besonders einfach. Ist nämlich (Bild 1) a ein 
Einheitsvektor in Richtung Q4 und m ein solcher 
in Richtung QM, so ist offenbar, wenn alle Summen 
vektoriell aufgefaßt werden, 

(AM)=(4AQ)-+ (QM) 
und daher 
—r AM)» — 4)+m(QM) 
zu, -4D* 2 nen il w? 


während mit a als Einheitsvektor La (Bild 1) 
offenbar: 





QAo!, _(QANe 
a = rg = j 


/ 





m 
In 
l 


— — (QM) wo? , - (Q4) 
AN, =b,- AM, =— m 2 + ., 


’ ‘ 


Ist ferner m ein weiterer Einheitsvektor _L m, so 
folgt wegen AN = AN, cot® und weil AN_ LAN, 
und aus AN, im Sinne von e gedreht werden muß 
mit Beachtung von (1) 


[ OM) w® Q, rn 
AN =| mM E Ba > < 


’ r Ja e 
—_ (QM) w»? (A) 
__ Mo’ _ QM 
ea a 
und daher wegen (2) 

; ww! — -- ( M) »° — —— > 
QN=QA+AN = m“ IK = — ım-QM- i 
ea d 


w.z.b. w.?) 


4. Allgemeine Folgerungen und spe- 
zielle Zerlegungen des Beschleunigungs- 


vektors. Aus den in Bild 1 angegebenen 
Koordinaten des Wendepoles W und Tangential- 
Bun ed 
poles T folgt PQ= === und somit wegen 
r J e? +1 2 
tgy= d [Gl. (2)] für die Koordinaten von @: 
5 Er. ed? 
Ig= PQsiniy=-, n s 
e? + d? 
(4). 
— e?d 
a Na PQ cos y=—— 
e? + d? 


Führt man (4) in die aus Bild 1 ersichtlichen 
Beziehungen, wobei in den Zählern die Koordinaten 
der Punkte Q,M,N erscheinen, nämlich in: 


Zy—79 %0-%u - Yy—YQ  Tum—2Q 
—— = m — Ä — == "—— 
ON Om ON Om 


ein und beachtet (2), so folgt nach einigen Kürzun- 
gen für die Koordinaten ry und ’yy des Punktes N 


e e z 
Ixv=e-—ı . YUyv=XT a Te 
N YJM d YN M d \ 


das nicht nur für die Zeichnungsgrößen, sondern 
ebenso auch für die natürlichen Abmessungen, wie 
man leicht einsieht, gilt, so daß in (5) die Querstriche 
sämtlich gleich fortgelassen wurden. 


Beschreibt M eine Gerade, so folgt dies auch 
für N, wie sich schon aus dem Begriffe der Dreh- 
streckung ergibt. Wandert insbesondere M auf 
der y-Achse, so bewegt sich N auf der z-Achse, 
denn aus (5) folgt für zy=0 sofort yy=0. Ist 
überdies „y=0, so folgt zy=e, d.h., wenn M in 
den Drehpol P fällt, so liegt N im Tangentialpol 7 
und wir haben die eingangs erwähnte W.Schell- 
sche Zerlegung der Beschleunigung vor uns. Ist 
aber neben ry=0 noch 4y =d, so folgt aus (5) 
neben 4y =0 sofort zy=(. d. h. die Koinzidenz 
von M mit dem Wendepole W bedingt die von N 
mit dem Drehpele P und wir haben die eingangs 
erwähnte M. Grüblersche Zerlegung erhalten. 


,) 


a d i 
Ist schließlich yy=— , 80 ist = woraus 


> 


der Satz folgt: 


Jede’Systembeschleunigung läßt sich 
in zwei Komponenten zerlegen, von 
denen die eine als relative Normalbe- 
schleunigung gegen den Mittelpunkt des 
Wendekreises, die andere als relative 
Tangentialbeschleunigung gegen den 
Mittelpunkt des Tangentialkreises auf- 
gefaßt werden kann. 


5) Ein analytischer Nachweis, der sich ebenfalls un- 
schwer führen läßt, möge hier wegen Raumersparnis 
unterdrückt werden. Ein anderer synthetischer Beweis 
rührt, wie ich einer frdl. schriftlichen Mitteilung ent- 
nehmen konnte, von Herrn Dr. R. Beyer her: Er wird in 
der Zeitschr. „Getriebetechnik‘ mitgeteilt werden. 
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Bild 2. Spezielle Zerlegung des Beschleunigungsvektors 
b4a:1./inn4aMo ttaN. ?inn4T+ttiak. 


In Bild 2 ist diese Zerlegung der Beschleunigung 
AA, = 4Qı + AR, durchgeführt: Q,M.ı || QM, er- 
gibt nm AM., während M„A,ı = AN. =t ano 
ist.und daher vom Mittelpunkt N, des Tangential- 
kreises unter dem <{® erscheint, wobei, tg9 = 
AR, . —. ; 

70 ist; außerdem muß AN,, 1 AN, sein, was als 
weitere Kontrolle dienen kann. 
man ferner leicht den Satz: 

Beschreibt M eine z- (y)-Parallele, so 
beschreibt N eine %- (r-)Parallele und 
umgekehrt. 

Beschreibt ferner M einen Kreis, so tut dies 
auch N, wobei die Mittelpunkte nach (5) trans- 








Aus (5) folgert 


= 
formiert und die Halbmesser im Verhältnis ; ver- 
( 


größert werden. Beschreibt insbesondere M den 
Wendekreis. so beschreibt N den Tangentialkreis. 
Beschreibt M den Rückkehrkreis, sp beschreibt N 
den Kreis mit dem Durchmesser TS, der durch 
Spiegelung des TangentialKreises an seiner Tangente 
durch T entsteht und der samt dem Rückkehrkreis 
in Bild 2 strichliert eingezeichnet wurde.‘ Liegt M 
insbesondere im Rückkehrpol R, so liegt N im 
Punkte S auf der Abszissenachse, somit in der 
Entfernung 2e vom Drehpol P, ein Ergebnis, das 
auch sofort aus (5) gefolgert werden kann. 

Liegt schließlich M im Tangentialpol 7, so liegt N 


i R e®? * i 
im Punkte K (e. 7); der im Schnitte der %-Par- 
C 


allelen durch 7 mit der Verbindungslinie PQ von 
Dreh- und Beschleunigungspol erhalten wird. 
Liegt jedoch N im Wendepol W,. so liegt M im 


(d? { ’ i 2 
Punkte L(- .d) ,„ der wieder im Schnitt der z- 


Parallelen durch W mit PQ gefunden wird. 

Alle diese Zerlegungen lassen sich leicht durch 
Sätze analog dem oben für die Zerlegung M,— N, 
ausgesprochenen formulieren, wobei man bloß die 
Bezugspunkte entsprechend einzusetzen hat. Natür- 
lich können derartige Zerlegungen in beliebiger 
Zahl angegeben werden. Die Zerlegung der Be- 
schleunigung 

b4 = AA, BR AT, +T,4;; (T, 4, . AK,) 
in naar und £4x istin Bild 2 eingezeichnet worden, 
wobei die in Bild 2 der Deutlichkeit halber nicht 
gezeichnete Gerade Q, T, || Q7 ist. 
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Bild 3. Zur Konstruktion der Beschleunigungskomponenten 


n4N und £ 4 M und deren Summe nn tt4M- 








Bild 4. Zur Lösung der Aufgabe: A zu suchen, wenu 
<I(mAam,stan=e, <Z[(b4,vA)=P 
vorgeschrieben sind und M auf k, N auf g liegen soll 


Zum Schlusse möge bemerkt werden, daß man 
auch die relative Tangentialbeschleunigung ?4m 
von 4 geren M und die relative Normalbeschleuni- 
gung n ıy von A gegen N leicht finden kann: Ihre 
Summe ist freilich nicht d 4, sondern b4—-dby-br 
wie man folgendermaßen leicht zeigen kann: Es 
bestehen die beiden Gleichungen: 


ba=bytnamttaım 
ba=bytnan titan 





2b, =butbntnanttan)tnanttan: 


Im Additionsergebnis ist der Klammerausdruck, 
wie bewiesen, aber selbst b ,, so daß 


ba— bu —bn=nanttan 


folgt, wobei natürlich sämtliche Additionen vek- 
toriell zu verstehen sind. Außerdem ist d5y — ba 
(daQ@M_L_QN) und By fällt in die Richtung NM. 
Es gelten somit die Zerlegungen: 


n +tın=b 
AM AN 4 | (6). 


NANT tim = ba a by _— dx j 
In Bild 3 ist n,x durch Ziehen von Parallelen 
zu QN bzw. MN aus n;g oder nyy direkt kon- 
struiert worden, während ?;y aus ?4n oder £4Q 
nach deren Drehung um 90° in AN bzw. AQ und 
nach Ziehen der Parallelen zu NM bzw. QM und 
Rückdrehung des Ergebnisses um 90° (Methode 
der gedrehten Beschleunigungen) gefunden wurde. 
Schließlich ist aus 5b ,, mittels der Methode der 
gedrehten Beschleunigungen dy und by gefunden 
worden, deren Subtraktion von db; in der Tat in 
den Endpunkt der Vektorsumme n4y-+!.m führt, 
wie es nach (6) sein muß. 
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6. Beispiel: Es sollen jene Systempunkte A 
gesucht werden, für die <(n4m,t4ıx)=«a und 
<(54,v4)=B vorgeschrieben ist, während M auf 
einem gegebenen Kreise Ak, N aber auf einer ge- 
gebenen Geraden g liegen soll. Man konstruiert 
zunächst aus %k (Bild 4) mittels der durch den 
Beschleunigungswinkel y gegebenen Drehstreckung 
den Kreis kA’, der g in 2 Punkten schneidet. Nach 
Wahl des einen derselben (N) findet man rück- 
wärts M auf k. Der Systempunkt 4 muß offenbar 
auf jenem Kreis k, durch M und N liegen, der den 


Peripheriewinkel = —+ a besitzt, und außerdem 
auf jenem Kreis k, durch P und Q, der den Peri- 
besitzt. Die Kreise A, 


und Ak, schneiden sich im allgemeinen in 2 Punkten, 
so daß günstigstenfalls 4 Lösungen der Aufgabe 
existieren ®). 


pheriewinkel y+#3 — = 


Prag. Karl Karas. 522 


6) Man vgl. hinsichtlich des Kreisesk3 L. Burmester: 
Lehrbuch der Kinematik, Leipzig 1888, S. 807 und Fig. 780 
(im zugehörigen Atlas, 3. Lieferung, Tafel 52). 


Zur Ingenleur-Axonometrie. Bei: der Her- 
stellung anschaulicher Skizzen technischer Gegen- 
stände benutzt der Ingenieur mit Vorteil die „In- 
genieur-Axonometrie“, d. h. eine senkrechte Pro 
jektion, bei der sich die Verkürzungszahlen e,, 
ey, 'e, der Koordinatenachsen z, y. 2 wie 1:2:2 
verhalten‘), Die Einfachheit dieser Darstellung 


beruht — vor allem beim Skizzieren von Dreh- 
körpern — darauf, daß sich für die Verkürzungs- 


zahlen ea, e3, ey der drei Koordinatenebenen 
atz.32__y,»r_.z einfache Werte ergeben’), 
nämlich 


10’ Br 


ea 3° 
(Das folgt sofort aus der Forderung 


ex:&y:e,=1:2:2 
und den Gleichungen 
e te, +e=2, tte=e,+e=e+e=1.) 


Für die Konstruktion dieses Achsenkreuzbildes 
benutze ich seit Jahren in meinen Vorlesunger: 
statt der üblichen Konstruktionen von Paster- 
nak oder von C. Volk?°) die folgende ein- 


1) Verkürzungszahl einer Geraden oder einer Ebene = 
Cosinns ihres Neigungswinkels gegen die Bildebene 

2) Blatt DIN5 des Deutschen Normenausschusses. 

3) ©. Volk: Die maschinentechnischen Bauformen und 
das Skizzieren in Perspektive. Berlin 1942. Pasternak, 
Ens. Math., XXIV, 1924 bis 1925. 


N 








Bild 1. Achsenbild der Ingenieuraxonometrie. 





fachere Überlegung. Es sei XYZ das (zunächst 
noch unbekannte) Spurendreieck, 0° der Höhen 
schnittpunkt, also das Bild des Koordinaten 
anfangspunktes OÖ, und N der lHlöhenspurpunkt im 
rechtwinkligen Dreieck XOY (Bild 1). Die Pro- 
jektion ist dimetrisch, also 0’Y —0'Z. Setzen 
wir ON —1, so ist im rechtwinkligen Dreieck 
NOZ, dem „Profildreieck“ der Z-Achse: Kathete 
ON =3-O’N =3, Hypotenuse NZ=9, also 02 
=0'Y=8. Nun folgt im Dreieck NZY, daß NY 
a 
=y63, also YZ=12=-,0'2 
Damit ergibt sich folgende Vorschrift: Bei der 
Ingenieur-Axomometrie. wählt man 
zunächst auf dem Bilde der Z-Achse 
eine beliebige Strecke O’Z und kon- 
struiert dann den Punkt Y so, daß 





0Y=0'Z und ZY= . :-0'Z 








Schließlich wird das X-Achsenbild 
lYZskizziert,also.als Halbierungs- 
liniedes stumpfen Winkels zwischea 
Y- und Z-Achsenbild. Die Seiten des 
Spurendreiecks werden natürlick nicht mitge- 
zeichnet. 

Bei dieser Gelegenheit sei darauf hingewiesen, 
daß man sich zum wirklichen Einzeichnen von 
Gegenständen in dieses Achsenkreuz am zweck- 
mäßigsten ein für allemal eine Profilfigur her 
stellt, die (ohne Benutzung eines Spurendreiecks!) 
so entsteht (Bild 2): 

1. Man zeichnet das Profildreieck NOZ der Z 
Achse mit O’N =1 (beliebig!), ’Z-8, NO=3. 
Es ist zugleich das Profildreieck der Y-Achse. 


Bild-Maßstäbe: wahrer Moßstab 
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Bild 2 Profilfigur der Ingenieuraxonometrie. 
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2. Aus:dem wahren Maßstab werden mit Hilfe 
dieses Profildreiecks die Bildmaßstäbe für die Y- und 
Z-Achse und für die Neigungslinien der f- und y- 


s . Zn l 

Koordinatenebenen | Verkürzungszahl ) gewon- 
2 / 

nen. (Die Konstruktionen sind in der Figur durch 

Pfeile angedeutet.) 

3. Aus dem Z-Achsen-Bildmaßstab wird durch 
Verkleinern auf die Hälfte der X-Achsen-Bildmaß- 
stab hergestellt und daraus das Profildreieck LOX 
der X-Achse (siehe Bild 1). 

4. Aus diesem ergibt sich der Bildmaßstab für 


die Neigungslinien der Koordinatenebene a | Verkür- 


9 
zungszahl x 25) 

Diese Profilfigur wird so benutzt: Für Strecken. 
die parallel sind zu den Koordinatenachsen oder 
parallel zu den Neigungslinien der Koordinatenebe- 
nen und deren Maßzahlen bekannt sind, werden die 
Bildmaßstäbe benutzt. Sind jedoch diese Strecken 
nicht durch ihre Maßzahlen, wohl aber in einer 
Grund- und Aufrißzeichnung gegeben, so greift 
man z. B. die Strecke r, die die Richtung der Z 
Achse habe, mit dem Stechzirkel aus der ge- 
gebenen Zeichnung ab, trägt sie von O aus auf 
OZ bis P ab, läßt die Zirkelspitze in P und ver- 
kürzt die Zirkelöffnung so, daß die zweite Spitze 
auf O0’ schleift, hat also die verkürzte Höhe im 
Zirkel. Ebenso wird z. B. durch Benutzung von 
( 
pr verkürzt, also die 
Nebenachse .der Ellipse gewonnen, die das Bild 
eines in a liegenden Kreises mit dem ‚Radius r ist. 
Durch Benutzung dieser (am besten vorgedruckten) 
Profilfigur werden also die sonst üblichen Um 
legungen der Koordinatenebenen und die Umrech- 
nung von Maßen vermieden. 

Braunschweig. Fritz Rehbock. 534 


OL ein Radius r auf > 


Netztafel-Nomogramme aus orthogo- 
nalen Kreisscharen. Die nomographische 
Darstellung einer Funktion F(a,db,c)=0 wird 
man, wenn durch Netztafeln, dann durch solche 
aus Geradenscharen zu erhalten trachten, schon 
aus Gründen einfacher Herstellung. Bequem her- 
zustellen sind aber auch Netztafeln mit einer, zwei 
oder drei Kreisscharen. Das Grundsätzliche über 
diese Tafelform sowie über ihre an Beispielen durch- 
geführte Herstellung findet man u. a. bei Werk- 
meister!). Wie aber die dort S. 33 gegebene 
allgemeine Gleichung zeigt, handelt es sich 
zumeist um verwickelter gebaute Funktionen 
F (a,b,c) — 0, die durch solche Tafeln zu erfassen 
sind. Bei den einfacheren Fällen ist oft zugleich 
eine Dä:stellung z. B. durch Geradentafeln mög- 
lich, denen dann der Vorzug zukommt. Bekannt- 
lich benützt man Netztafeln aus 3 Kurvenscharei 
so, daß man den Schnitt einer bezifferten Kurve 
der ersten Schar mit einer solchen der zweiten 
Schar sucht und dann die Bezifferung derjenigen 
Kurve der dritten Schar abliest. die durch den 
Schnitt der beiden ersten geht. Das setzt einen 
scharfen Schnitt der beiden ersten Kurven voraus. 
Dies ist bei Netztafeln mit Kreisscharen nicht 
immer und i. a. auch nicht auf allen Teilen der- 
selben Tafel gleichmäßig der Fall, ein weiterer 
Grund, die Anwendbarkeit der Nomogramme aus 
Kreisscharen herabzusetzen. 

Der letztgenannte Übelstand wird aber grund- 
sätzlich beseitigt, wenn man Tafeln benutzt, bei 
denen zwei der Kreisscharen zueinander ortho- 


1, P. Werkmeister: Das Entwerfen von graphischen 
Rechentafeln. Berlin 1923. 


gonal sind./ Gibt man ihnen die Bezifferung der 
beiden unabhängigen Veränderlichen, so ist durch 
ihren immer einwandfreien Schnitt die Kurve der 
dritten Schar genau zu legen und abzulesen. Wie 
Werkmeister a. a. O. S. $0ff. ausführt, be- 
nutzt man zum Zeichnen der drei Scharen ein 
Hilfssystem, in dem die Gleichungen der 3 Scharen 
ausgedrückt werden. Die Elimination der Hilfs- 
veränderlichen x und y ergibt dann F (a,b,c) — 0, 
eben die durch das Nomogramm dargestellte 
Funktion. 

Im folgenden werden zwei orthogonale Kreis- 
scharen (in einfacher und symmetrischer Lage zum 
Hilfssystem) mit. einer dritten, ebenfalls einfachen 
und einfach gelegenen Geradenschar bzw. Kreis- 
schar kombiniert und daraus nach obigem Ver 
fahren die jeweils dargestellte Funktion, F (a,b, c) 

0 ermittelt. Statt der Parameterwerte . der 
Scharen, d. h. der dargestellten Veränderlichen 
selbst, wird man in der Endgleichung natürlich 
Funktionen dieser Veränderlichen f (a), g (b), h (c) 
setzen können. 


A. Die in Bild 1 gezeichneten orthogonalen 
Kreisscharen haben die Gleichungen 
AI: 2 +? —2ar — #?=0 
AU: 2+ 7 —2by+p=I, 
wobei » die im Bild angegebene Konstante ist. 
1. Wählt man als 3. Schar die Ursprungsstrahlen 
mit der Gleichung =cx, so ergibt die Differenz 


der Gleichungen AI und AII in Verbindung mit 
dieser Gleichung die Werte 


p? 


ne SE. SE. 
a y 


—cb—a' 
Diese werden in AI (oder AII) eingesetzt, und es 
entsteht nach einigen Umformungen (mit dem Ziel 
der möglichst symmetrischen Form) die durch das 
Netz dargestellte F(a,b,c)=0 in der Gestalt 

a? —= c? (b? — 9?) — p?. 


2. Die 3. Schar sei jetzt die der Parallelen zu 
einer der Achsen. Derselbe Lösungsweg gibt jetzt 
im Falle: 


yac bzw. g=cC 
die Werte 


Das Einsetzen in AII bzw. AI liefert die darge- 
stellte Funktion 
(+ 0) (pP — 2bc)+c?(a +b)—=0 
für die Schar „= c bzw. 
(Pb) (pP +2ac)+c (a +b)=0 
für die Schar 2 = c. 


3. Die 3. Schar bestehe aus Kreisen um den Ur- 
sprung mit der Gleichung 2 +9? —c=0. Das 
Einsetzen von c? für (2?+y?) in AI und AI 
ergibt die Werte 


c? + p® 
u 
und deren Einsetzen in AI und AII führt nach 


einigen Umformungen zur Gleichung der darge- 
stellten Funktion 


c? Eee p? 
& 52 und y= 


Ma.E 2 2 x 
(c* Ha + 2.) +2 pw e(l, — .,) —4°=0. 
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Bild 1 und 2. 


Bild 3 (unten). Die Funktion 
b_i+e 


a 1—e’ 








Diese Funktion, ebenso wie die unter 2., wird in 
der Praxis selten vorkommen, so daß sich in 
diesem Betracht die Herstellung solcher Nomo- 
gramme wohl nicht lohnen würde, trotz ihrer er- 
wähnten Genauigkeit und ihrer verhältnismäßig 
einfachen Herstellbarkeit. 


(Kreismittelpunkte Halbmesser 
Sobar I: zu =, ya=0 ?—=p!+ a? 
Schar I: z2,=>0, Yamb ”=b?—-p. 


Die. Halbmesser sind nach diesen Gleichungen 
auch leicht geometrisch zu finden, wie Bild 1 an- 
deutet.) 


B. Wenn man in den vorstehenden Endgleichun- 
gen p=0 setzt, entsteht aus dem ÖOrthogonalnetz 
von’ Bild 1 dasjenige des Bildes 2, und aus den 
Gleichungen unter A 1, 2 und 3 werden bezüglich 
die folgenden Gleichungen unter B 1, 2 und 3, die 
man nach Jem genannten Verfahren auch unmittel- 
bar aus den Gleichungen der Scharen von Bild 2: 


BI: #2-+y —2ar=0 
BI: +7 —2by=0 


und einer dritten Schar herleiten kann. Die sich 
ergebenden einfachen Funktionen seien nun be- 
trachtet. 


1. Die 3. Schar besteht aus den Ursprungs- 
strahlen y= cz. Dies Nomogramm bringt Werk- 
meister a. a. 0. S. 53. Es stellt einfach die 
Funktion a=b:c dar. Sind a und 5b die unab- 
hängigen, c die abhängige Veränderliche, so wird 
der Ursprung mit dem Schnitt je eines a-Kreises 
und eines b-Kreises verbunden (wirklich oder bloß 
durch angelegtes Lineal) und der zugehörige 
Wert c, sei es auf einer mit c bezifferten Geraden 
parallel der y-Achse, sei es auf einem mit c be- 
zifferten Kreis um den Ursprung, abgelesen. Statt 
der Veränderlichen a,d,c können Funktionen der- 
selben, f(a), g(b), h(c), gesetzt werden, so daß 
entsteht: ’@ —=g(b)-h(c). Diese Funktion wird 
ebensogut dargestellt durch eine Strahlentafel; 
zwei Scharen /(a) und g(b) sind die Achsen- 
parallelen, die Werte A (e) sind die Steigungen der 
Ursprungsstrahlen. 


2. Die 3. Schar besteht aus den Parallelen zu 
den Achsen, und zwar folgt für die Schar 


y=c bzw. #>° 
als dargestellte Funktion: 


sn 2ab , b 
— at D? ZW, 


2ab? 
Tarn 


Ablesebeispiel Punkt P: 


q=1 db=3 c=üß. 














oder umgeformt: 


1 1 2 1 1 2 
= —, bzw. - ee, 
a? ” bs Dec a? + b? ac 
Setzt man noch im ursprünglichen Ansatz oder in 
der Endgleichung 2c für c, so kommt 


we Ka ER TR 
a? + b? be bzw. a? + b ac 


als verhältnismäßig einfache, durch das Nomo- 
gramm dargestellte Funktion. Der große Vorzug 
(neben der Ablesegenauigkeit) ist die leichte Her- 
stellbarkeit. Es sind 


die Kreismittelpunkte die Halbmesser 


von Sebar I: z=e, y=d Pe 
von Schar II: z=0, y=b tb. 


3. Die 3. Schar seien Kreise um den Ursprung: 
2”+y?=c’. Die unmittelbare Rechnung oder 
das Nullsetzen von p» in der Gleichung unter AS 
liefert: 

Es 1 4 
a? 


DD 0 
oder, wenn man 2c statt .c setzt: 


1 1 1 
„r»”7 ; 


c? 
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Dies Nomogramm übertrifft in der Einfachheit der 
Herstellung und Anwendung noch das vorige. Zu 
beachten ist besonders, daß die Mittelpunkte bzw. 
Halbmesser der beiden ÖOrthogonalscharen linear 
auf der Achsen aufgetragen werden, nicht qua 
dratisch oder reziprok quadratisch, wie es bei 
anders gebauten Nomogrammen für dieselbe 
Funktion der Fall sein würde. 


4. Es sei noch das Nomogramm mit denselben 
Örthogonalscharen besprochen, wenn die 3. Schar 
aus Parallelen zur 135°-Linie oder zur 45 °-Linie 
besteht. Deren Gleichungen sind: 

y+tı=c bzw. y—ı=c. 
Aus‘der Differenz BI—BIl folgt mit diesen 
Gleichungen: 


dass b c Y — 
a—+b' ; a-+b 
bzw. 
__ be BE... 
Tr ai Te 


Durch Einsetzen in BI und nach einigen Um- 
formungen folgt, mit 2c statt c, 


1 1 Eur. -) 
ar bc (3 a a, 
bzw. 
1 1 1/3 1 
at bc (5 ai ) 
als Gleichung der dargestellten Funktion. 


C. Eine ähnlich einfache Lage haben die beiden 
Orthogonalscharen, wenn ihre Zentralen um 45° 
geneigt werden (Bild 3). Ihre Gleichungen sind 
dann 


CI: 2 +9 —2a(r+y)=0 
CN =+y—2b(r —J)=0. 
a und 5b haben jetzt die im Bild ersichtliche Be 


deutung. Als dritte Schar kann man wieder die 
bisherigen Scharen benutzen. 


1. Die dritte Schar besteht aus den Ursprungs- 
strahlen y=czr. Die Differenz CI— CII gibt in 
Verbindung damit die dargestellte Funktion: 


1+c 
G-— 


1—c' 


b—o 


.b+a 


c oder b= 


te gleich A(c), so hat man wieder 
die Funktion unter Bi: b=-@c mit der Strahlen- 
tafel als einfacherer Darstellungsform. Behält 


Setzt man 


man aber c bei oder setzt für c die Funktion f (c), 


so hat das Kreisnomogramm bezüglich Herstellung 
entschieden einen Vorrang wor der Strahlentafel. 


Denn bei dieser ist der Ausdruck ei eine Stei- 


gung. Man könnte ihre Werte auf einer Ge- 
raden parallel zur y-Achse auftragen, müßte dann 
aber diese Gerade projektiv nach c teilen; nicht 
iinear. Beim Kreisnomogramm wird dieselbe Ge- 
rade Hinear beziffert. Ein ausgeführtes Beispiel 
zeigt Bild 3. 


2. Als dritte Sehar benutzt man die Achsen- 
parallelen 


y=c bzw. 


Die Rechnung gibt 


b+a } 
z=cC - )ZW, = 
b a 3 


b ee 
b-a 


und weiter, mit 2c statt c, als dargestellte Funk 
tion: 


1 1 Br 1 R r 

a? r be c\a b 
bzw. 

1 1 wu) 

at u” (7 7%)’ 


also im wesentlichen dieselbe Funktion wie im 
Fall B4. 


3. Die 3. Schar seien die Kreise um den Ur- 
sprung: 2?+y?=c’. Ersetzen von (z?+?) 
durch ce? in CI und CII gibt 2 Gleichungen in z 
und % mit den Lösungen 

ce a-+b ce a—b 


x = — - —— 
5 und Y 


deren Einsetzen in CI ergibt: 


1 1 8 
-— 0. 


a? b? c? 


PETER 
art 


naturgemäß dieselbe Funktion wie unter B3. 


4. Die 3. Schar bestehe aus Parallelen zur 135 °- 
bzw. zur 45°-Linie, mit den Gleichungen 


y+r=c bzw. y—-ı=c. 
Die Herleitung liefert 
© c 
=> t 4 = —— { — 
z 3, Pre, Y 5, ® a) 
bzw 
c c 
ı=— ;.(e+rb), y-5,(@-8. 


Das Einsetzen in CI bzw. CI ergibt, mit 4c 
statt c: 


u =% SR 
tr bo ac B. er” be 


im wesentlichen dieselbe Funktion wie die unter 


Die Betrachtung zeigt. daß Netztafel - Nomo- 
gramme, bei denen zwei Scharen aus orthogonalen 
Kreisen in einfacher symmetrischer Lage zum zu- 
grundeliegenden Koordinatensystem bestehen, 
leicht herstellbar sind. Verwendet man sie so, 
daß die beiden unabhängigen Veränderlichen durch 
die beiden Kreisscharen dargestellt werden, so ist 
die Genauigkeit des Nomegramms groß, denn die 
zugehörige dritte. Veränderliche steht an einer 
Linie, die durch den scharfen Schnittpunkt der 
beiden Orthogonalkreise gut festgelegt ist. Die 
durch ‚solche Nomogramme dargestellten Funk- 
tionen sind, wenigstens bei den Kreisscharen- der 
Bilder 2 und 3, verhältnismäßig einfach gebaut, 
was für praktische Anwendungen von Bedeutung 
ist. 


Flensburg. Christian Ahrens. 5% 
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BUCHBESPRECHUNGEN 


Dr. LUDWIG PRANDTL, Prof. a. d. Univ. Göttin- 
gen, Führer durch die Strömungslehre, 
zugleich dritte Auflage des Abrisses der Strö- 
mungslehre. XI+383 S. m. 314 Abb. Braun- 
schweig 1942, Verlag Friedr. Vieweg & Sohn. 
Preis geb. 12 M. 


Der Prandtische „Abriß der Strömungslehre“, 
der nunmehr in 3. Auflage vorliegt, verdankt 
seine Sonderstellung unter der zeitgenössischen 
Lehrbuchliteratur der Tatsache, daß-er trotz des 
Verzichtes auf längere mathematische Deduktionen 
nicht bei den elementaren Gegenständen stehen 
bleibt, sondern den Leser in klarer und auschau- 
licher Weise mit den Fragestellungen und Ergeb- 
nissen der weitverzweigten modernen Strö- 
mungsforschung vertraut macht. Die erste Auf- 
lage des Abrisses erschien 1931, die zweite kam 
1935 heraus und war im wesentlichen ein unver- 
änderter Abdruck der ersten. Dagegen ist in der 
jetzt erschienenen dritten Auflage der Umfang 
des Buches stark angewachsen, von 223 auf 382 
Seiten, so daß die Bezeichnung des Werkes als 
Abriß nicht mehr angemessen erschien. 


Einmal sind in den schon bisher vorhandenen 
ersten vier Abschnitten bedeutsame Moderni- 
sierungen und Ergänzungen angebracht worden. 
Um den umfassenden Standpunkt des Werkes an 
einem ersten Beispiel zu belegen, sei etwa er- 
wähnt, daß im Anschluß an die technischen Pro- 
peller nunmehr auch die im Tierreich auftretenden 
Propulsionsarten behandelt werden. Vor allem ist 
aber ein fünfter groß angelegter Abschnitt über 
Randgebiete der Strömungslehre ganz neu in das 
Buch aufgenommen worden. In loser Anordnung 
behandelt dieser Abschnitt Forschungsgebiete, die 
etwas isoliert von der üblichen Strömungslehre 
stehen. Die Fachwelt ist dem Verfasser nicht nur 
für die mühselige Sammlung und Überarbeitung 
dieses weitschichtigen Materials zu Dank ver- 
pflichtet, sondern auch für die erstmalige Mit- 
teilung zahlreicher eigener Überlegungen, die neue 
Aufklärung bringen. Eine ganze Anzahl dieser 
Gegenstände hat nahe Beziehungen zu Frage- 
stellungen der Geophysik, so der Sandtransport 
durch den Wind und die Geschiebebewegung in 
Flüssen. Speziell mit der dynamischen Meteoro- 
logie sind eng verknüpft die Kapitel über Strö- 
mungen in geschichteten schweren Flüssigkeiten 
und über die Strömungen in rotierenden Bezugs- 
systemen, die aber auch für Strömungsmaschinen 
von erheblicher Bedeutung sind. Technische Son- 
dergebiete betreffen die Ausführung über Kavi- 
tation, Wasserschlag und Gleitflächen. Ferner 
wird das verwickelte Gebiet der aufgezwungenen 
und der natürlichen Strömungen bei ungleich- 
mäßiger Temperaturverteilung behandelt, also der 
Fragenkreis, der mit dem Wärmeübergang zusam- 
menhängt. Als. „natürlich“ wird dabei eine erst 
durch die Temperaturverteilung selbst erzeugte 
Strömung bezeichnet. Mit diesen für die Technik 
wie für die Geophysik in gleicher Weise wichtigen 
Darlegungen schließt das imposante Werk, das 
wie kein zweites das universale Eingreifen der 
modernen Strömungslehre in naturwissenschaft- 
liche und technische Fragestellungen zeigt t). 


Dresden. W. Tollmıen. 547 


1) Die im Juli 1944 erschienene, vom stehen gebliebenen 
Satz gedruckte vierte Auflage, weist gegenüber der dritten 
nur geringe Änderungen auf. Es handelt sich dabei 
um kleinere textliche Verbesserungen und Ergänzungen, 
sowie einige Hinweise auf inzwischen erschienene Ver- 
öffentlichungen. W. 


P. RAETHJEN, Einführungindie Phy- 
sikderAtmosphäre,Bd.Il, Meteorolo- 
gische Aerodynamik. (Hamburger Mathe- 
matische Einzelschriften, 33. Heft/1942.) VIII+ 
257 S. m..57 Abb. Leipzig und Berlin 1942, Verlag 
PR, G. Teubner. Preis geb. 10 M. 


Auch in diesem zweiten Band seiner Einführung 
in die Physik der Atmosphäre, der die meteoro- 
logische Aerodynamik behandelt» hat der Ver- 
fasser mit großem pädagogischen Geschick deü 
sehwierigen Stoff gemeistert. Vor. die mathema- 
tische Deduktion, die mit einfachen Beweismitteln 
die allgemeinen Gesetzmäßigkeiten liefert, wird 
meist die Beschreibung schlichter Experimente ge- 
stellt, durch welehe die Fragestellung sinnfällig 
gemacht wird. Übungsaufgaben dienen dem Leser 
zur Überprüfung des Verständnisses. 

Der Stoff ist in zwei Kapitel geteilt, in klein- 
räumige und großräumige Aerodynamik. Bei der 
kleinräumigen Aerodynamik wird die Corioliskraft 
nicht berücksichtigt. In diesem Kapitel werden 
also aus der gewöhnlichen Aerodynamik diejeni- 
gen Kenntnisse vermittelt;die auch für die meteorolo- 
gischen Anwendungen grundlegend sind. Gleich- 
zeitig wird aufgewiesen, welche meteorologischen 
Vorgänge sich noch ohne Berücksichtigung der 
Erdrotation behandeln lassen. Die Thermodynamik 
spielt bei den Ausführungen über Cumulusaufwind 
und Thermikaustausch sowie über Feuchtlabilitäts- 
energie hinein. Die der dynamischen Meteorologie 
eigentümlichen Züge kommen allerdings erst recht 
in dem Kapitel über großräumige Aerodynamik 
zur Geltung. Der Aufbau dieses reichhaltigen 
Kapitels sei durch die Paragraphenüberschriften 
gekennzeichnet: Rotationsmoment und Trägheits- 
kreis; Coriolisbeschleunigung und reibungslose 
Strömungsgleichungen; die bodennahe Reibungs- 
schicht; Strömungsgleichgewicht der freien Ät- 
mosphäre; Wetterluftmassen, Frontalzonen, Grenz- 
flächen; Stabilität und Labilität des geostrophi- 
schen Windfeldes;' Zirkulationsbeschleunigung. 

Es ist gewiß, daß dieses Buch durch seine an- 
genehme Darstellungsart und die ständige Bezug- 
nahme auf die meteorologische Praxis diesem 
schwierigen, aber zukunftsreichen, Zweig der Strö- 
mungslehre neue Freunde erwerben wird. 


Dresden. W. Tollmien. 548 


Dr. sc. techn. ERWIN HOECK, Druckver- 
luste in Druckleitungen großer 
Kraftwerke. Bericht über -die Arbeiten der 
Druckverlustkommission des Schweizerischen In- 
genieur- und Architektenvereins. (Mitteilungen aus 
der Versuchsanstalt für Wasserbau an der Eidg. 
Techn. Hochschule in Zürich Nr. 3.) 75 S. m. 
45 Abb. Zürich 1943, Verlag AG. Gebr. Leemann 
& Co. Preis Fr. 10,— +4% Umsatzsteuer und 
Porto. 

In der vorliegenden Mitteilung wird ausführlich 
über ausgedehnte Messungen des Druckverlustes 
in den Druckleitungen zahlreicher schweizerischer 
Kraftwerke berichtet. Eine ausführliche Diskussion 
wird in der Richtung geführt, dem Ingenieur ver- 
besserte Mittel zur Eiulsung neuer Anlagen an 
die Hand zu geben. Sie erfolgt an Hand der 
Widerstandsformeln für rauhe Rohrleitungen von 
Prandtl-Nikuradse, Che6zy, Strick- 
ler, Ludin. Es zeigt sich, daß moderne, voll- 
ständig geschweißte Rohrleitungen (mit Ausnahme 
anz neuer spiegelglatter) bei den praktisch in 

tracht kommenden großen Re ein quadratisches 
Widerstandsgesetz Jiefern, daher sowohl nach Ni- 
kuradse, wie auch nach Ch&zy oder Strick- 
ler berechnet werden können, wenn man die äqui- 





Bemz» zum 












Z. angew. Math, Mech. 
Band 24 Nr. 2/1944 


Buchbesprechungen 





9: 




















valente Sandrauhigkeit bzw. die Widerstands- 
koeffizienten aus den berichteten Messungen ent- 
nimmt, wofür eine geeignete Tabelle und ein ent- 
sprechendes Kurvenblatt gegeben werden. Da- 
gegen bieten Leitungen mit genieteten Quernähten 
und geschweißten oder genieteten Längsnähten 
mehr Schwierigkeiten. Hier kann nach den Mes- 
sungsergebnissen die Formel von Ludin v=k 
.rU;65.. 70,58) am besten verwendet werden. Alles in 
allem genommen muß die sehr gründlich durch- 
geführte Arbeit als ein wichtiger Beitrag zum 
Problem des technischen Rohres großer Abmes 
sungen gewertet werden, wenn auch in Anbetracht 
der vielgestaltigen Möglichkeiten die sichere Klassi- 
fizierung nach wie vor -weiter Schwierigkeiten 
bieten wird. 

Leipzig. L. 


Schiller. 5832 


Dr.-Ing. habil. KARL KLOTTER, Prof. a. d. 
Techn. Hochschule Berlin, Messung mecha- 
nischer Schwingungen (Dynamik der 
Schwingungsmeßgeräte). IX +151 S. m. 
80 Abb. Berlin 1943, Springer-Verlag. Preis brosch. 
8,40 M. 

Die Bücher über Meßgeräte befassen sich vor- 
wiegend mit beschreibenden Darstellungen ihrer 
Wirkungsweise und ihrer, Handhabung, während 
die theoretischen Unterlagen nur in Einzelfällen 
wiedergegeben werden. Das. vorliegende Buch be- 
handelt‘ in systematischem Aufbau die Dynamik 
der Geräte zur Messung mechanischer Schwin- 
gungen. Nach einem kurzen Abriß der Schwin- 
gungslehre, bei der von der anschaulichen vekto- 
riellen Darstellung weitgehend Gebrauch gemacht 
wird, werden die für Meßgeräte wichtigen Ver- 
größerungsfunktionen, Phasenverschiebungswinkel 
und Phasenverschiebungszeiten dargestellt. Mit 
diesen Grundlagen werden die Geräte zur Kraft- 
und Bewegungsmessung, die auch die Geschwin- 
digkeits-, Beschleunigungs- und Rückmessung mit 
umfaßt, hinsichtlich ihrer Abstimmung und zweck- 
mäßigen Dämpfung erläutert. Im einzelnen werden 
bei der Bewegungsmessung die Anordnungen mit 
und ohne Festpunkt, sowie mit Feder- und mit 
Reibungsfesselung erörtert. Im letzten Abschnitt 
des Buches werden die Verfahren zur Entzerrung 
bei nicht-periodischen Einwirkungen beschrieben. 

Das Buch schließt eine merkliche Lücke im 
Fachschrifttum und bildet für jeden mit Aufgaben 
des Messens betrauten Wissenschaftler ein not- 
wendiges Rüstzeug. Auch der Schwingungstheo- 
retiker wird in dem mit großer Sorgfalt ausgear- 
beiteten Buch, dessen geschlossene Darstellungs- 
weise allerdings das auszugsweise Nachlesen ein- 


zelner Abschnitte erschwert. eine Reihe neuer An-. 


regungen finden. 


Stuttgart. P. Riekert. 529 


‚Dr. se. techn. EL-SAYED EL-SCHASLY, Bie 
gungsspannungen und Stabkräfte in 
Schwedlerkuppeln nach Theorie und 
Modellversuch, miteinerAÄnwendung 
aus dem Flugzeugbau. (Mitteilungen aus 
dem Institut für Baustatik an der Eidg. Techn. 
Hochschule Zürich, herausgegeben von Prof. Dr. 
M. Ritter und Prof. Dr. F. Stüssi, Nr. 12.), 185 S. 
m. 76 Fig. u. 3 Tafeln. Zürich u. Leipzig 1943, Ver. 
lag A. G. Gebr. Leemann & Co. Preis brosch. 9 M. 


Die Fehier, die sich durch die Annahme von 
reibungslosen Gelenken an Stelle der biegesteifen 
Anschlüsse der Stäbe räumlicher Fachwerke im 
allgemeinen und bei der Schwedlerkuppel im be- 
sonderen ergeben, sind durch mehrere Beiträge 
zur Berechnung räumlicher Fachwerke mit zykli- 
scher Symmetrie bekannt. Man weiß, daß sich der 


Spannungszustand der idealen Schwedlerkuppel 
mit reibungslosert Gelenken wesentlich von dem 
Kräftebild der ausgeführten Stahlkonstruktion 
unterscheidet. Der Verfasser hat sich in seiner 
Arbeit die dankenswerte Aufgabe gestellt, die ein- 
zelnen Annahmen zur statischen Untersuchung 
derartiger Tragwerke, die Schwedlerkuppel mit 
gelenkigen Stabanschlüssen, die Schwedlerkuppel 
mit biegesteifen Sparren, die Schwedlerkuppel mit 
biegesteifen Ringen, die Bedeutung des biege- 
steifen Laternenringes für den Spannungszustand 
und die Schwedlerkuppel, mit biegesteifen Spar- 
ren und Ringen, insbesondere mit biegesteifen 
Sparren und biegesteifem Schlußring zu unter- 
suchen, an Hand eines Zahlenbeispiels ‘quantitativ 
zu klären und damit die Bedeutung der Annahmen 
zur Vereinfachung der Rechnung zu bewerten. Die 
Brauchbarkeit der Rechnung ist darauf durch die 
Messung von Spannungen und Knotenverschie- 
bungen an einem Zeluloidmodell geprüft und 
durch die statische Untersuchung der Rumpfnase 
eines mehrmotorigen Flugzeugs mit gutem Erfolg 
gezeigt worden. Die Arbeit zeichnet sich durch 
Sorgfalt und Sachkenntnis aus. Sie ist als Disser- 
tation der Eidgenössischen Technischen Hoch- 
schule Zürich erschienen und kann als Vorbild 
für ähnliche Untersuchungen aus dem Fachgebiet 
der Baüingenieurwissenschaften dienen, 


Dresden. K. Beyer. 543 
Dr. ABD-EL-AZIZ EL-AROUSY, Studien 
über das elastische Verhalten von 


Brückengewölben einschließlich des 
Zusammenwärkens mit dem Aufbau. 
(Mitteilungen aus dem Institut für Baustatik an der 
Eidg. Techn. Hochschule Zürich, heratisgegeben von 
Prof: Dr. M. Ritter u. Prof. Dr. F. Stüssi, Nr. 13.) 
132 S. mit 23 Abb. Zürich u. Leipzig 1942, Verlag 
A.G. Gebr. Leemann & Co.- Preis brosch. 7,20 M. 


Das elastische Zusammenwirken eines Bogen- 
trägers mit dem aufgesetzten oder angehängten 
Fahrbahnrost ist in der deutschen technischen 
Literatur mehrfach und zuerst wohl von 
Schachenmayer im Jahre 1908 in einer Arbeit 
über mehrfache Brückengewölbe untersucht wor- 
den. Die Verbindung der beiden Baukörper besitzt 
ohne Zweifel für die strenge Beschreibung des 
elastischen Zusammenhangs und damit für die 
zuverlässige Beurteilung der Festigkeit von Bogen- 
träger und Fahrbahnrost Bedeutung. Die Unter- 
suchung wird indessen in der Reget unterlassen, 
da der‘ Verschiebungszustand eines eingespannten 
Bogenträgers unter der Wirkung der Nutzlast 
klein ist und vernachlässigt wird, also die Festig- 
keit des Fahrbahnrostes für die Nutzlasten unter 
der Annahme eines starren Bogenträgers mit un- 
verschieblichen Anschlüssen der Pfosten unter- 
sucht wird. Sie führt außerdem auf die Bearbei- 
tung eines hochgradig statisch unbestimmten Trag- 
werks und verlangt die Auflösung eines weit- 
gehend zusammenhängenden fehlerempfindlichen 
Systems von linearen Gleichungen. Die Zusatz- 
spannungen aus dem Zusammenwirken von Bogen 
und Fahrbahnrost gelten dann als Nebenspannun- 
gen. Diese werden von dem Verfasser nach ver- 
schiedenen Ansätzen qualitativ und quantitativ an 
Hand zahlreicher Tabellen und Beispiele sorgfältig 
untersucht, so daß der Leser ein gutes Bild über 
den Einfluß des elastischen Zusammenhangs-auf 
den Spannungszustand des Tragwerks gewinnt. 
Darin liegt der Wert der Veröffentlichung. Das 
Studium der Arbeit bietet jedenfalls dem Leser, 
insbesondere dem. Entwurfsverfasser von Bogen- 
brücken, wertvolle Anregungen und kann aus 
diesem Grunde bestens empfohlen werden. 

K. Beyer. 544 


Dresden. 
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ZUSCHRIFTEN AN DEN HERAUSGEBER 


F. Sehubert, Zur Theorie des stationären 
Verdichtungsstoßes. Z. angew. Math. Mech. Bd. 23 
(1943), S. 129 bis 138. 

Die in obiger Arbeit hergeleiteten Formeln, die 
die Zustandsgrößen hinter einem schrägen Ver- 
dichtungsstoß als Funktion von M, siny darstellen, 
wobei M, die Machsche Zahl vor dem Stoß und y 
den 'Neigungswinkel des Stoßes gegen die An- 
strömrichtung bezeichnet, sind bereits in einem 
französischen Lehrbuch enthalten. Im Jahre 1921 
ist von L. Jakob ein Buch „La Rösistance -de l’Air 
et l’Experience — Less Cönsequences“, Tome I im 
Verlag Doin, Paris, erschienen, in dem auf S. 250 
die Formeln (45) bis (51) zusammengestellt sind, 
die mit denen von F. Schubert genau überein- 
stimmen. Lediglich in Formel (50) ist im Nenner 
ein Summand vergessen worden, wie auch Tesson 
erwähnt hat, der in seiner Arbeit „Contribution 
ä l’etude des ondes de derangement“ in den 
Comptes Rendus, Bd. 217, 1943, S. 208 auf obiges 
Buch hingewiesen hat. 

Braunschweig. E. Groth. 

Erwiderung. Das von Herrn Groth erwähnte 
Buch von L. Jakob ist mir nicht bekannt. Dem 
Titel nach dürfte es sich weniger um ein Lehrbuch 
als um eine Monographie handeln. Da ich .das 
Buch von Jakob bisher nicht einsehen konnte, 
war es mir noch nicht möglich, die Herleitungen 
von L. Jakob mit den meinigen zu vergleichen. 
Ein wesentlicher Gesichtspunkt meiner Arbeit war 
aber gerade die möglichst einfache und durchsich- 
tige Beweisführung. Daß manche der Formeln als 
solche schon anderweit mitgeteilt waren, hatte ich, 
soweit mir dies noch vor der Veröffentlichung be- 
kannt wurde, schon angemerkt (vgl. Fußnote 3 
meiner Arbeit). Jedenfalls hatte sich die vorteil- 
hafte Parameterdarstellung der Verdichtungsstoß- 
theorie nicht durchgesetzt. Meine Arbeit war aus 
der Praxis der Berechnungen kompressibler Strö- 
mungen: erwachsen und sollte dieser Praxis wieder 
dienen, wofür nicht zuletzt die ausführlichen 
Zahlentafeln und ee Diagramme Zeugnis 
ablegen. F. Schubert. 558 


H. Karl,‚Biegung gekrümmter, dünnwandiger 
Rohre. Z. angew. Math. Mech. Bd. 23 (1943), S. 331 
bis 345. 

Im Rahmen meiner Mitarbeit in einem Sonder- 


ausschuB des Kriegsmarine - Normenausschusses 
habe ich mich vor einigen Jahren mit theoretischen 
Untersuchungen über die Beanspruchung und die 
Formänderung von glatten Rohrkrümmern infolge 


von Biegemomenten befaßt, die in bzw. senkrecht 
zu der Ebene des Krümmers wirken. Ich bin dabei 
abweichend von den Kärmänschen Ableitungen 
nicht von der Formänderung, sondern von den 
Gleichgewichtsbedingungen am Krümmerelement 
ausgegangen und habe für die Spannungen in der 
Umfangsrichtung des Krümmers (Längsspannun- 
gen) genau die gleichen trigonometrischen Reihen- 
entwicklungen gewählt, wie sie auch in der Arbeit 
von Herrn Dr.-Ing. Karl benützt werden. Die 
Koeffizienten dieser Reihen habe ich dann in be- 
kannter Weise aus den Bedingungen für das Mini- 
mum der Formänderungsarbeit ermittelt; im Hin- 
blick auf die für praktische Anwendungen er- 
forderliche Genauigkeit habe ich mich mit der 
zweiten Näherung begnügt. In meinen Ableitungen 
ist die Verschiedenheit der Krümmungsradien ver- 
nachlässigt; den Einfluß der Spannungen in der 
Umfangsrichtung des Rohrquerschnitts (Tangen- 
tialspannungen) sowie der Querkontraktion habe 
ich- nur bei Ermittlung der Werte für die erste 
Näherung berücksichtigt. Die Ergebnisse ' meiner 
Untersuchungen waren entsprechend den grund- 
sätzlich gleichen Ansätzen völlig die gleichen wie 
diejenigen, die Herr Dr.-Ing. Karl in den Ab- 
schnitten BIb, BII und CI seines Aufsatzes 
erhalten hat. 

Eine Veröffentlichung meiner Ausarbeitungen ist 
bisher noch nicht erfolgt; die Ergebnisse derselben 
sind jedoch bereits längere Zeit vor Einreichung 
der Dissertation des Herrn Dr.-Ing. Karl (im De- 
zember 1941 bzw. Januar 1942) mehreren Mit- 
gliedern des oben genannten Sonderausschusses 
bekanntgegeben worden. 


Darmstadt. R. Wälde. 

Erwiderung. Die Anwendung des Castigliano- 
schen Minimalprinzips zur Bestimmung der Bie- 
gung von Rohrkrümmern durch Biegemomente in 
der Krümmerebene und senkrecht zu dieser war 
durch Arbeiten von Lorenz und Weber bekannt. 
Im Gegensatz zur nichtveröffentlichten Arbeit von 
Herrn Prof. Wälde wird in meinem Aufsatz nicht 
nur die zweite Näherung gebracht, sondern die 
allgemeine Lösung für beide Minimalprinzipien und 
damit eine Eingrenzung und Fehlerbestimmung. 
Außerdem wird noch für die allgemeine Lösung 
der Einfluß der Querzusammenziehung und des 
Krümmungsverhältnisses untersucht und die Über- 
einstimmung der Koeffizienten für beide Arten von 
Biegemomenten für das Krümmerelement nachge- 
wiesen. 


Dresden. H. Karl. 56% 


NACHRICHTEN 


Persönliches. 

Der Dozent für reine und angewandte Mathematik 
an der Techn. Hochschule Berlin, Dr. phil. habil. 
Günther Schulz, wurde zum außerplanmäßigen 
Professor ernannt. 

Dr.-Ing. habil. Karl Marguerre wurde unter 
Belassung am Institut für Festigkeit der Deutschen 
Versuchsanstalt für Luftfahrt zum a. o. Prof. im 
Reichsdienst ernannt. 

Der Dozent an der Techn. Hochschule Karlsruhe, 
Dr. habil. L. Collatz, wurde am 31. 12. 1943 
als Nachfolger des verstorbenen Prof. Dr.G. Prange 
zum o. Prof. der Mathematik an der Techn. Hoch- 
schule Hannover ernannt. 

Als Nachfolger von Prof. M. Lagally wurde der 
em. ord. Prof. der Bergakademie Freiberg Dr. 


Fr. A. Willers zum ord. Prof. der angewandten 
Mathematik an der Techn. Hochschule Dresden er- 
nannt. 

Der em. o. Prof. der angewandten Mathematik 
an der Techn. Hochschule Dresden, Dr. M. La- 
gally, wurde zum Mitglied der 'mathematisch- 
physikalischen Klasse der sächsischen Akademie 
der Wissenschaften gewählt. 

Die vor 50 Jahren vom VDI gestiftete Graß- 
hof-Denkmünze, die als erster ©. von Bach 
erhielt, wurde nach mehrjähriger Pause in diesem 
Jahr wieder verliehen, und zwar an Prof. Dr. 
A. Thum. 

Bei einem Luftangriff kam der Dozent an der 
Techn. Hochschule Braunschweig, Dr.-Ing. habil. 
Walter Hantzsche, ums Leben. 
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